BROADBAND  PUMPING  EFFECTS  ON  THE 
DIODE  PUMPED  ALKALI  LASER 


THESIS 


Paul  Nathaniel  Jones,  Civilian 
AFIT/GAP/ENP/1 1-M04 

DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright-Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  refleet  the  offieial 
policy  or  position  of  the  United  States  Air  Force.  The  Department  of  Defense,  or  the  United 
States  Government.  This  material  is  declared  the  work  of  the  U.S.  Government  and  is  not 
subject  to  copyright  protection  in  the  United  States. 


AFIT/GAP/ENP/1 1-M04 


BROADBAND  PUMPING  EFFECTS  ON  THE  DIODE  PUMPED 

AEKAEI  EASER 

THESIS 


Presented  to  the  Faeulty 
Department  of  Engineering  Physies 
Graduate  Sehool  of  Engineering  and  Management 
Air  Eoree  Institute  of  Teehnology 
Air  Einiversity 

Air  Edueation  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Seienee  in  Applied  Physies 


Paul  Nathaniel  Jones,  BS 
Staff  Seientist,  AFIT 


Mareh  2011 


APPROVED  FOR  PUBEIC  REEEASE;  DISTRIBUTION  UNEIMITED. 


AFIT/GAP/ENP/1 1-M04 


BROADBAND  PUMPING  EFFECTS  ON  THE  DIODE  PUMPED  ALKALI  LASER 


Paul  Nathaniel  Jones,  B.S. 
Civilian 


Approved; 


_ //signed// _  _ 

Glen  P.  Perram,  PhD  (Chairman)  Date 


_ //signed// _  _ 

Lt.  Col  Jeremy  C.  Holtgrave  (Member)  Date 


_ //signed// _ 

Maj.  Clifford  Sulham  (Member) 


Date 


AFIT/GAP/ENP/1 1-M04 


Abstract 

This  research  seeks  to  gain  greater  insight  on  the  mechanics  of  The  Diode  Pumped 
Alkali  Laser  through  analytic  modeling  techniques.  This  work  is  an  extension  to  a 
previous  model  developed  by  Dr.  Gordon  Hager,  and  focuses  on  the  addition  of  pump- 
beam  bandwidth.  Specifically,  it  seeks  to  determine  the  effect  that  broadband  pumping 
has  on  laser  performance.  The  model  incorporates  all  the  fundamental  parameters  within 
the  laser  system,  including  alkali  concentrations,  collision  partner  concentrations,  pump 
bandwidth,  length  and  temperature  of  gain  medium,  transmission,  and  reflectivity. 
Baseline  operating  conditions  set  Rubidium  (Rb)  concentrations  ranging  from  10^^  -  10*"^ 
atoms/cm  ,  corresponding  to  operating  temperatures  ranging  from  50  -  150  C.  Ethane  or 
Methane  concentrations  are  varied  corresponding  to  partial  pressures  from  100  -  500 
Torr.  The  system  is  evaluated  for  incident  beam  intensity  ranging  from  0  -  30000 
W/cm  ,  for  both  lasing  and  non-lasing  system  analysis.  Output  laser  beam  intensities 
scale  well  with  input  beam  intensity  and  the  model  predicts  optical  to  optical  efficiencies 
of  over  70%. 
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BROADBAND  PUMPING  EFFECTS  ON  THE  DIODE  PUMPED  AEKALI  EASER 


I.  Introduction 

Background 

For  decades  the  United  States  Military  has  been  researching  and  developing  high 
energy  lasers  (HEE)  for  improved  tactical  advantages  and  anti-ballistic  missile  defense. 
These  weapons  offer  numerous  advantages  over  conventional  systems,  and  as  Secretary 
of  the  Air  Force  Widnall  recognized  in  1997:  “It  isn’t  very  often  an  innovation  comes 
along  that  revolutionizes  our  operational  concepts,  tactics  and  strategies.  You  can 
probably  name  them  on  one  hand  -  the  atomic  bomb,  the  satellite,  the  jet  engine,  stealth, 
and  the  microchip.  It’s  possible  the  airborne  laser  is  in  this  league.”  With  pinpoint 
accuracy,  HEE  are  capable  of  reducing  collateral  damage  while  achieving  target 
destruction.  Their  light  speed  trajectory  enables  them  to  engage  targets  immediately  after 
detection  from  a  relatively  safe  distance.  And  finally,  they  are  inexpensive  to  operate  and 
can  be  tailored  for  other  purposes,  including  non-lethal  destruction,  and  disruption  [8]. 

The  currently  mounted  HEE  on  the  Airborne  Easer  (ABE)  has  proven  effective  yet 
possesses  a  few  undesirable  characteristics.  In  February  2010,  The  ABE  flew  over  the 
coast  of  California  and  successfully  targeted,  tracked,  and  destroyed  a  ballistic  missile  in 
its  boost  phase.  This  historical  feat  marks  a  breakthrough  in  operational  status  and 
revolutionizes  warfare.  In  fact,  the  Missile  Defense  Agency  officially  recognized  this 
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fact  by  awarding  its  Technology  Pioneer  Award  to  three  Boeing  Airborne  Laser  Testbed 
engineers  and  three  of  their  government  and  industry  teammates.  Yet  while  this 
achievement  is  groundbreaking  and  a  warfare  game-changer,  eriticism  over  the  installed 
systems  logistical  feasibility,  size,  and  ammunitions  reservoir  is  justified.  This  system  is 
a  chemieally  driven  laser  whieh  utilizes  vast  quantities  of  hazardous  chemieals  ineluding 
chlorine,  iodine,  hydrogen  peroxide  and  potassium  hydroxide,  none  of  whieh  are  ideal  for 
storage  in  a  war  zone.  After  exhausting  its  6-10  shot  magazine  depth,  it  would  take  two 
C-17  transport  planes  to  refuel,  and  due  to  atmospheric  turbulence  the  ABL  has  a  limited 
range  of  approximately  600  km  for  liquid  fueled  ICBM  and  300  km  for  solid  fueled 
ICBM.  In  addition,  it  takes  the  entire  hull  of  a  Boeing  747  to  ereate  weapons  grade 
power.  However,  it  is  still  the  only  system  eapable  of  such  high  power  on  target.  With 
these  factors  in  mind  a  pursuit  for  a  lightweight,  eompact,  low  hazard  risk,  and  high 
magazine  depth  system  has  been  underway. 

The  Diode  Pumped  Alkali  Laser 

According  to  experts  in  the  field,  the  system  which  holds  the  most  promise  for 
aehieving  these  standards  is  the  Diode  Pumped  Alkali  Laser  (DPAL).  It  was  first 
demonstrated  in  the  early  90 ’s  and  patented  by  Krupke  in  2003  [5].  These  systems 
utilize  the  cost  and  energy-effectiveness  of  diode  bars,  the  thermal  management 
advantages  of  gas  phase  lasers,  the  simplieity  of  alkali  earth  metals,  and  are  essentially  an 
efficient  narrow-banding,  laser  photon  engine  [2]. 

An  alkali,  typieally  Rubidium  (Rb)  or  Cesium  (Cs),  is  heated  to  its  gas  phase  and 
subject  to  diode-driven  laser  light  in  resonance  with  its  second  lowest  excited  state  (^P3/2). 
A  spin-orbit  relaxer  induces  a  rapid  Boltzmann  Equilibrium  between  the  upper  pump 
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state  to  the  upper  lasing  state,  or  the  lowest  lying  exeited  state  (  P1/2).  Spontaneous 
emission  from  the  P1/2  state  to  ground  induces  amplified  stimulated  emission  leading  up 
to  lasing.  The  energy  level  diagram  below  shows  the  optical-atomic  and  kinetic 
processes  involved. 


Energy  Level 

’3/2^— 


7\ 


Spin-Orbit  Relaxation/Excitation 
I_ i 


Pump  Absorption  /  D2 
(780  nm) 


S.E.  /  Laser  Beam  /  D1 
(795  nm) 


Figure  1.  DPAL  Three  Level  Diagram  (Rb) 


The  DPAL  mass  produces  single  mode,  narrow-banded,  near  infrared  (IR)  laser 
photons  in  a  cyclic,  three  step  process.  As  shown  above,  broadband,  multimode  diode 
laser  light  centered  at  the  D2  transition  is  incident  on  a  vaporized  alkali  and  buffer  gas 
mixture.  This  light  excites  the  alkali  into  its  second  excited  state.  With  the  proper  buffer 
gas  to  alkali  concentrations,  rapid  Boltzmann  Equilibrium  occurs  through  collision 
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processes  causing  relaxation  from  the  upper  P3/2,  to  the  lower  P1/2  state  [2],  While  the 
energy  split  between  the  upper  pump  and  upper  lasing  level  is  less  than  1%  of  the  pump 
energy,  at  projected  operating  temperatures  this  equilibrium  still  places  over  45%  of  the 
alkali  atoms  in  the  upper  lasing  state,  resulting  in  a  high  population  inversion.  Once  the 
alkali  is  in  this  state,  (the  upper  lasing  level),  it  slowly  (with  a  radiative  lifetime  of  27.7 
ns  [5]),  yet  eventually  begins  to  fluoresce  back  to  ground.  This  causes  those  slowly 
trickling  spontaneous  emission  photons  to  see  an  inverted  landscape  upon  release.  Once 
ejected,  they  pass  through  the  gain  medium  and  by  stimulated  emission  “grab”  photons 
from  other  excited  Di  atoms.  As  they  grab  these  photons  the  projecting  wave  gains  in 
amplitude  through  amplified  spontaneous  emission  (ASE)  until  reaching  the  end  of  the 
gain  cell. 

Given  a  laser  cavity  and  an  emission  trajectory  within  a  mode  of  the  respective 
cavity,  and  assuming  our  gain  exceeds  threshold,  this  amplified  wave  will  oscillate  back 
and  forth  through  the  gain  medium  exponentially  increasing  in  amplitude  after  each  pass. 
Eventually  and  after  multiple  passes,  the  stimulated  emission  rate  caused  by  the  high 
intensity  oscillating  laser  beam  overcomes  the  pump  and  spin-orbit  relaxer  rates,  leading 
to  a  saturated  Di  to  ground  transition  and  the  maximum  intra-cavity  laser  intensity  [2]. 
Analytical  Modeling 

In  order  to  accurately  model  these  events  and  predict  the  output  beam  strength  of  the 
DEAL,  a  robust  and  quasi-analytic  model  has  been  developed  by  Dr.  Gordon  Hager  [2]. 
This  model  accounts  for  all  the  substantial  interactions  and  effects  associated  with  DEAL 
operation,  including  input  beam  intensity,  alkali  concentration,  buffer  gas  concentration, 
mirror  reflectivity,  cavity  transmitivity  and  gain  medium  length.  It  integrates 
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fundamental  laser  physies  with  experimental  results  speeifie  to  the  DPAL  to  produee 
feasible  and  aeeurate  results. 

It  has  been  advaneed  in  order  to  take  into  effeet  the  frequeney  dependenee  of  the 
pump  beams  interaetion  with  the  system.  It  ereates  an  alkali  eross-seetion  bandwidth 
funetion  that  is  dependent  on  temperature  and  pressure,  and  also  adds  a  frequeney 
dependenee  to  the  pump  beam  intensity.  It  shows  that  a  elose  speetral  overlap  between 
these  two  funetions  is  paramount  for  optimal  effieieney.  Finally,  it  is  also  a  benefieial 
tool  for  analysis  of  gain  medium  kineties,  parametrie  studies,  and  helps  to  improve  the 
overall  intuition  of  the  staff,  in  hopes  that  a  more  elear  and  fundamental  understanding 
will  lead  to  an  effieient,  powerful,  and  preeisely  eontrolled  system. 
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II,  Literature  Review 


In  this  chapter,  a  review  of  literature  pertaining  to  the  physics  involved  with  the 
DPAL  is  summarized.  Also  summarized  is  a  review  of  publications  on  advancing  DPAL 
technology  and  a  few  small  references  to  laser  engineering.  After  establishing  the 
theoretical  foundations  of  the  DPAL  and  its  current  state  of  affairs  an  introduction  to  the 
analytic  model  is  made  in  the  following  chapter  by  a  review  of  two  of  its  author’s 
publications. 

Alkali  Earth  Metals 

Alkali  earth  metals  (alkalis)  have  been  extensively  researched  and  are  an  ideal 
candidate  for  an  optically  pumped  gain  medium.  First,  they  consist  of  the  elements 
possessing  a  single  valence  electron,  making  their  electron  energy  levels  relatively  easy 
to  assign  and  well  documented.  The  National  Institute  of  Science  and  Technology  has 
publicized  a  table  consisting  of  the  more  significant  spectroscopic  constants  for  each 
element,  and  an  excerpt  from  the  table  for  Rb  is  shown  below: 


Table  1.  Rb  Spectroscopic  Data  [5] 


Wavelength  (nm) 

AtiS* 

Configurations 

Terms 

Ji  Jk 

gi  gk 

740.81731 

-i-p^5p  -  4p*7s 

2P-»  -  2S 

3  2  -  12 

4-2 

761.8933 

4p*5p  -  4p*5d 

2P*  -  2D 

12-32 

2-4 

775.76507 

4p®5p  -  4p*5d 

2P*  -  2D 

3  2  -  5  2 

4-6 

775.94363 

4p®5p  -  4p®5d 

2P*-2D 

3  2  -  3  2 

4-4 

780.0268 

3.81x10* 

-l-p®5s  -  4p®5p 

2S  -2P» 

12  -  32 

2-4 

792.526 

4p®4d  -  4p®8f 

2D  -2F* 

5  2  -  7.2 

6-8 

792.5537 

4p*4d  -  4p*8f 

2D  -2F* 

3  2  -  52 

4-6 

794.7603 

3.61x10’ 

4p*5s  -  4p^5p 

2S  -2P* 

12  -  12 

2-2 

827.141 

4p*4d  -  4p“7f 

2D  -2F* 

5  2  -  7,2 

6-8 

827.1707 

4p^4d  -  4p^7f 

2D  -2F* 

32-5  2 

4-6 

886.8512 

4p*4d  -  4p*6f 

2D  -2F* 

52-72 

00 

886.8852 

4pNd  -  4p*6f 

2D  -2F* 

3  2  -  5  2 

4-6 
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A  tabulation  of  each  electronic  transition  in  Rb  corresponding  to  photonic 
wavelengths  from  740-880  nm  along  with  the  key  spectroscopic  constants  of  Einstein  A 
coefficients  (column  2),  initial  and  final  electronic  configurations  (column  3),  and  the 
degeneracy’s  of  the  initial  and  final  states  (column  5).  Key  transitions  in  the  DEAL  (Di 
at  795  nm  and  D2  at  780  nm)  have  well  known  wavelengths,  Einstein  A-coefficients  , 
degeneracy  ratios,  and  spin  states.  Eor  a  more  complete  listing  of  data  for  Rb  see 
Appendix  A.  The  absorption  cross-section  is  dependent  on  the  Einstein  A-coefficient  and 
is  given  by: 

-  -  (1) 

Where  g2  is  the  degeneracy  of  the  upper  transition  state  and  gi  of  the  lower.  is  the 
Einstein  A  coefficient,  is  the  transition  wavelength,  and  is  the  normalized,  unitless, 
lineshape  function  (typically  a  Gaussian,  Eorentzian,  or  Voigt  profile).  It  is  apparent  that 
absorption  cross-sections  scale  linearly  with  Einstein  A  coefficients  and  it  should  also  be 
noted  that  alkalis  possess  some  of  the  largest  coefficients  currently  known,  at  linecenter 
the  Di  transition  absorption  cross-section  has  a  published  value  of  4.92x10'  cm  [2,5]. 
Given  such  high  cross-sections,  alkalis  are  easily  and  efficiently  energized  through  a 
radiant  source,  or  in  the  case  of  the  DEAL,  through  diode  laser  pumping. 

Another  characteristic  amongst  the  alkalis  that  proves  them  an  effective  laser  gain 
medium  is  their  extremely  high  quantum  efficiency  (QE).  In  a  laser,  QE  refers  to  the 
ratio  between  the  excitation  transition  and  the  lasing  transition.  Mathematically  stated: 

-  (2) 
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Where  is  the  QE, 


and 


,  are  the  pump  and  laser  transition  energies 


respeetively.  It  is  the  upper  limit  for  the  performanee  of  any  laser.  It  eannot  be  designed 
and  is  solely  dependent  on  energy  levels  of  the  gain  medium.  However,  lasers  ean  be 
designed  around  QE,  and  the  DEAL,  having  only  a  three  step  proeess  and  merely  a  fine 
structure  split  between  its  pump  and  upper  lasing  state,  has  remarkable  QE.  A  tabulation 
of  each  alkali’s  QE  is  shown  here: 


Table  2.  Alkali  Quantum  Efficiencies  [5] 


Alkali 

AEpump  (sV) 

AEiaser  (^V) 

QE 

Av  (GHz) 

H 

1.8887008 

1.8886887 

0.9999936 

2.9257441 

Ei 

1.84786 

1.847818 

0.9999773 

10.155475 

Na 

2.10442907 

2.10229704 

0.9989869 

515.51853 

K 

1.61711291 

1.60995774 

0.9955753 

1730.0989 

Rb 

1.589049 

1.5595909 

0.9814618 

7122.8813 

Cs 

1.45520606 

1.38648645 

0.9527767 

16616.198 

Er 

1.7263556 

1.5172452 

0.8788718 

50562.275 

Each  alkali  is  listed  with  its  photonic  pump  and  laser  energy,  the  QE,  and  the 
difference  in  frequency  (Av)  between  the  pump  and  laser  transition.  Such  high  QE  are 
ideal  for  optimal  performance  but  a  large  enough  frequency  gap  must  be  maintained  in 
order  to  avoid  spectral  overlap  between  the  absorption  and  lasing  transition.  As  will  be 
shown,  alkalis  under  current  experimental  conditions  have  absorption  cross-sections  with 
a  full  width  half  max  (EWHM)  on  the  order  of  100s  of  MHz,  leaving  plenty  of  room  (for 
most  of  the  alkalis)  in  frequency  space  to  avoid  significant  overlap.  However,  it  is  this 
authors  opinion  that  Rb  or  Cesium  (Cs)  are  the  most  ideal  candidates.  They  maintain 
very  high  QE  (  >  95%),  and  have  thousands  of  gigahertz  in  spacing  between  their 


transitions,  boosting  confidence  that,  even  at  high  pump  intensities  and  alkali 
concentrations,  any  lineshape  overlap  is  negligible. 

While  alkalis  possess  extremely  high  absorption  cross-sections,  and  have  remarkable 
QE,  as  mentioned  before,  matching  the  spectral  bandwidth  of  the  pump  source  to  the  D2 
absorption  profile  is  paramount  for  optimal  efficiency  [9].  In  order  to  accurately  model 
this  profile  any  broadening  or  shifting  of  its  hyperfine  structure  must  be  well  accounted 
for.  Figure  2  below  diagrams  the  hyperfine  structure  on  the  D2  transition  for  Rb; 
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Figure  2.  Rb  D2  Hyperfine  Structure  [13] 
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As  seen,  there  are  four  hyperfine  components  in  the  5  P3/2  state  and  two  components 
in  the  5  S1/2.  However,  the  selection  rules  of  AF  =  0,  ±1  allows  for  only  a  total  of  four 
transitions  within  the  D2  line.  Recently  published  data  by  Pitz  et  al  has  observed  and 
accurately  modeled  the  Doppler  and  pressure  broadening  rates  as  well  as  any  pressure 
shifts  associated  with  each  of  these  transitions  for  Cs.  It  gives  the  net  absorption  cross- 
section  as: 

-  -  (3) 

Where  gj  2J+1  is  fine  structure  degeneracy, 

is  linecenter  frequency  of  the  F’^F”  hyperfine  component, 

2  is  c/y,  A2iis  the  spontaneous  emission  rate, 

Sf’f"  is  the  hyperfine  line  strength  for  the  component, 

Gv(v,vf’  ^f")  is  the  Voigt  lineshape  centered  at  vf- ^f", 
fiso  is  relative  natural  abundance  of  Rb,  and 
fF"  is  the  statistical  distribution  of  population  among  F”  states. 

This  model  has  zero  fit  parameters  and  agrees  with  high  precision  experimental  data. 
However,  at  high  temperatures  and  pressures  the  cross-section  blends  into  a  single, 
pressure  broadened  line.  These  rates  are  determined  for  Cs,  a  very  similar  element  to  Rb. 
Therefore,  the  broadening  rates  determined  by  this  fit  are  used  to  approximate  a  single 
Lorentzian  for  the  D2  transition  absorption  cross-section.  At  sub-atmospheric  pressures 
however,  a  peak  cross  section  value  error  of  13%  has  been  reported.  This  error  decreases 
with  pressure,  and  the  DPAL  will  potentially  operate  above  a  single  atm  in  order  to 
maintain  high  spin-orbit  relaxation.  Such  accuracy  in  the  absorption  lineshape  will  allow 
for  a  Gaussian  diode  beam  to  be  narrow-banded  for  maximum  absorption  [9]. 
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By  narrow-banding  the  pump  beam  in  such  a  way  and  maintaining  proper  buffer  gas 
to  alkali  concentrations,  and  assuming  no  lasing,  the  initial  pump  rate  and  spin-orbit 
relaxation  rate  far  exceed  the  spontaneous  or  stimulated  emission  rate  on  the  lasing 
transition.  This  leads  to  a  saturated  D2  transition  and  the  upper  two  levels  at  their 
Boltzmann  distribution.  In  this  instance  the  kinetics  of  the  system  can  be  modeled  by  the 
following  equations: 

-  -  —  (4) 

(5) 

(6) 

Where  is  the  total  alkali  concentration,  are  the  populations  in  the  ground, 

upper  laser,  and  upper  pump  state  respectively,  are  the  degeneracies  associated 

with  those  states,  AE  is  the  energy  split  between  the  upper  two  levels,  k  is  the  Boltzmann 
Constant,  and  T  is  temperature.  In  this  instance  two  limits  are  being  achieved 
simultaneously,  the  “bleached”  pump  transition  limit,  and  the  quasi  two-level  (Q2L) 
limit.  The  bleached  limit  corresponds  to  a  fully  saturated  pump  transition  and  the  Q2L 
refers  to  the  dynamics  when  the  spin-orbit  relaxation  rate  far  exceeds  the  spontaneous  or 
stimulated  emissions  rate  on  the  lasing  transition.  Solving  this  system  of  equations  gives 
the  pre-lasing  populations  of  the  three  states  involved,  and  as  mentioned  before,  at  room 
temperature  places  over  45%  of  the  population  in  the  upper  laser  level,  while  leaving 
only  17%  in  the  ground.  Such  a  high  inversion  allows  for  gain  above  threshold  to  be 
achieved  relatively  easily.  Figure  3  on  the  following  page  shows  the  inversion  ratio 
between  the  P1/2  and  ground  level  as  a  function  of  cell  temperature: 
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Figure  3.  Initial  Population  Inversion  vs.  Alkali  Temperature 


As  shown,  alkalis  (specifically  Rb),  are  capable  of  high  inversions  through  a  combination 
of  three  factors:  favorable  degeneracy  ratios  between  states,  rapid  interaction  with  buffer 
gases,  and  a  significant  Boltzmann  differential. 

In  order  to  ensure  this  inversion  is  maintained  and  accurately  predict  gain  medium 
kinetics  precise  knowledge  of  the  dynamics  between  the  alkali  and  buffer  gas  is  requisite. 
This  dynamic  enters  the  laser  rate  equations  in  the  form  of  the  spin-orbit  mixing  rate 

coefficient,  k32  with  units  of  cm  /s.  The  product  of  k32  and  the  buffer  gas  concentration 

-1 

(E),  in  cm'  ,  yields  the  number  of  collision  induced  energy  transfers  per  atom  per  second. 
Its  role  is  shown  explicitly  in  equation  7,  which  is  a  portion  of  the  lasing  mode  rate 
equation  for  the  population  density  of  the  upper  lasing  state: 

—  “  (7) 


12 


As  shown,  the  spin-orbit  mixing  rate  coeffieient  defines  how  quickly  the  upper  two  levels 
of  the  system  equilibrate.  The  factor  of  two  in  equation  7  is  due  to  the  2:1  degeneracy 
ratio  between  the  hyperfme  states  of  5  P3/2  and  5  P1/2  in  Rb,  respectively.  Theoretically 
the  rate  constant  is  defined  as: 

-  (8) 

Where  Vavg  is  the  average  relative  velocity  of  the  collision  pair,  p  is  the  reduced  mass  of 
the  collision  pair,  (p  =  [(MRbMMe)/(  MRb+MMe)],  and  a  is  the  spin-orbit  relaxation  cross- 
section.  This  rate  has  a  published  value  of  5.2x10'^°  cmVs,  and  at  700  torr  of  buffer  gas, 
cycles  each  atom  through  the  lasing  process  many  millions  of  times  per  second  [16]. 

This  allows  for  relatively  low  alkali  concentrations  to  create  high  power  intracavity  beam 
intensities,  potentially  minimizing  thermal  management  issues. 

This  spin-orbit  kinetic  process  cannot  be  fully  analyzed  without  an  accurate  account 
of  the  alkali  concentration.  At  projected  operating  temperatures  alkalis  are  in  their  liquid 
phase,  and  do  not  reach  their  boiling  point  until  1032  K.  However,  in  liquid  phase  there 
is  still  a  small  percentage  of  atoms  in  their  gas  phase  due  to  the  slight  probability  of  Rb 
atoms  separating  from  the  lattice  by  means  of  constructively  interfering  phonons, 
localized  “hot  spots”  within  the  liquid,  or  any  other  thermally  induced  energy  influx.  A 
compilation  of  Rb  data  by  Steck  gives  a  rough  guide  to  this  percentage  in  terms  of  vapor 
pressure  in  Torr  and  is  given  by  equation  9  below: 

-  (9) 

where  Py  is  the  alkali  vapor  pressure  in  Torr  [13].  A  graph  of  this  pressure  near  projected 
operating  temperatures  is  shown  here: 
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Figure  4.  Rb  Vapor  Pressure  vs.  Temperature  [13] 
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A  final  prerequisite  for  the  feasibility  of  the  DPAL  to  produee  weapons  grade 
intensity  and  beam  quality  on  target  is  high  atmospherie  transmittivity.  Atmospherie 
attenuation  is  devastating  to  the  advantages  of  laser  quality  radiation,  and  ean  greatly 
reduce  the  effective  range  of  the  weapon.  Absorption  bands  of  the  atmospheres 
constituent  particles  lead  to  diminishing  intensity,  dispersion,  and  diffraction  causing  a 
relatively  weak  and  dispersed  radiative  strike.  However,  at  the  operating  wavelength  of 
the  Rb  based  DPAL  (795  nm),  there  is  a  transmission  window  with  a  transmittivity 
coefficient  of  over  80%.  That  figure  is  the  value  for  the  transmission  through  the  entire 
atmosphere,  holding  promise  that  even  higher  transmission  will  occur  for  the  projected 
operating  range  of  the  ABL  (~500  km).  Figure  5  shows  the  atmospheric  percent 
transmission  for  a  broad  range  of  wavelengths  given  a  Zenith  angle  of  20  degrees; 


Figure  5.  Atmospheric  Transmittivity 
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Alkali  earth  metals  have  been  ehosen  as  the  gain  medium  in  the  DPAL  for  a  wide 
variety  of  reasons.  First,  there  is  a  strong  foundation  of  knowledge  regarding  their 
fundamental  subatomic  structure.  Second,  they  possess  remarkable  QE,  high  absorption 
probabilities,  and  are  capable  of  creating  large  population  inversions.  Also,  each  atom  is 
capable  of  producing  millions  of  laser  photons  per  second  through  a  high  Boltzmann 
differential  and  accessible  spin-orbit  relaxation.  They  are  easily  vaporized  at  modest 
temperatures  allowing  for  gas  phase  thermal  management.  They  produce  laser  photons  at 
frequencies  within  atmospheric  windows,  leading  to  high  intensity  weapons  grade  power 
on  target.  Finally,  they  are  cheap,  abundant  and  relatively  easy  to  use,  (albeit  reactive  to 
atmospheric  exposure).  They  hold  promise  for  high  power,  compact,  efficient,  and  cost 
effective  systems. 

Diode  Pumped  Alkali  Laser  Technology 

The  DPAL  has  generated  much  attention  from  many  institutions  including  The  Air 
Force  Institute  of  Technology,  Lawrence  Livermore  National  Laboratory,  University  of 
New  Mexico,  Air  Force  Research  Laboratory  (AFRL)  Directed  Energy  Directorate, 
Krupke  Lasers  and  others.  With  so  many  entities  in  pursuit  of  realizing  the  potential  of 
this  system  much  progress  has  been  made.  It  should  be  noted  however  that  the  DPAL 
technology  is  still  in  its  early  stages  and  there  is  still  much  design,  engineering, 
experimentation  and  ingenuity  required  before  this  dream  becomes  a  reality.  Presented 
here  is  a  review  of  publications  pertaining  to  technological  advancements  of  the  DPAL. 

Although  gas  phase  lasers  are  inherently  less  effected  by  thermal  heating,  issues  still 
arise  and  a  gain  medium  flow  system  has  been  successfully  designed,  developed  and 
demonstrated  by  AFRL.  This  system  flows  the  alkali-buffer  gas  gain  medium  through 
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the  resonator  of  the  DPAL  and  carries  away  heat.  By  doing  so,  a  fresh  set  of  atoms  is 
ready  for  the  absorption-relaxation  phases  of  the  3 -step  cycle  and  atoms  lost  to  effects 
such  as  multi-photon  absorption,  collisional  excitation  into  higher  lying  states,  ionization 
or  other  processes  outside  the  3  level  model  are  flown  out  of  the  cavity  [11],  This  is  a 
significant  first  step  towards  achieving  excellent  beam  quality  at  high  laser  powers. 

A  high  efficiency  Rb  DPAL  has  been  demonstrated  with  a  linear  optical  slope 
efficiency  of  69%.  In  this  setup,  a  static  alkali/ethane/helium  lasing  vapor  mix  is  pumped 
with  an  approximately  1  MHz  bandwidth  diode  laser  achieving  a  maximum  output  power 
of  490  mW.  This  was  achieved  with  400/50  Torr  of  ethane/helium  and  a  20%  reflective 
output  coupler.  It  also  achieved  an  optical-optical  efficiency  of  3 1 .5%  with  a  large 
threshold  of  -677  mW.  Linear  scaling  on  the  efficiency  was  maintained  showing  the 
spin-orbit  relaxation  rate  was  high  enough  to  maintain  a  quasi  two  level  system,  resulting 
in  a  proportional  increase  to  laser  intensity  as  input  beam  intensity  is  increased  [8]. 

A  second  demonstration  of  a  DPAL  device  by  Sulham  et  al  focusing  on  operating 
high  above  threshold  and  achieved  an  output  power  of  over  32  times  threshold  while 
maintaining  a  linear  slope  efficiency  between  39%  and  49%.  This  places  only  3%  of  the 
efficiency  loss  due  to  threshold  values.  It  also  suggests  that  no  second  order  kinetic  or 
optical  processes  are  involved  in  this  experimental  setup,  and  with  an  output  intensity  of 
over20kW/cm  demonstrates  relatively  high  output  powers  [13].  However,  the  pump 
beam  is  pulsed  with  a  100  ns  pulse  width  and  may  not  accurately  represent  continuous 
wave  operation,  since  thermal  effects  would  not  be  readily  apparent.  It  is  possible  that 
second  order  or  nonlinear  effects  may  not  set  in  until  a  greater  time  span.  Figure  6  at  the 
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top  of  the  page  1 8  shows  the  baste  setup  for  the  statie  eell  DPAL  experiments  mentioned 
above. 


An  alternative  approaeh  to  thermal  management  has  reeently  been  proposed  by 
Zweibaek  et  al  that  utilizes  a  high  aspeet  ratio,  liquid  eooled  eopper  waveguide  whieh 
serves  to  both  guide  the  pump  beam  and  more  quiekly  dissipate  heat.  In  traditional 
DPAL  designs,  the  volume  of  the  gain  cell  far  exceeds  the  actual  gain  volume  and  heat 
removal  is  only  accomplished  through  conduction  with  the  cells  walls.  The  large  spacing 
between  the  active  region  and  the  cell  wall  causes  low  thermal  conductivity  resulting  in  a 
thermally  constricted  upper  limit  on  laser  power.  In  this  design,  a  copper  core  runs 
through  the  gain  volume  dissipating  heat.  It  is  also  capable  of  heating  depending  on  the 
state  of  lasing.  Figure  7  at  the  top  of  the  follow  page  shows  the  lasing  head  and  gives 
output  power  as  a  function  of  time  [17]; 
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Figure  7.  Copper  Waveguide  Model  and  Results  [17] 


As  seen  on  the  left,  an  offset  coolant  channel  penetrates  the  gain  medium  rapidly 
dissipating  or  adding  heat  to  the  gain  medium.  The  plot  to  the  right,  as  interpreted  by  its 
author,  is  best  described  by  thermally  induced  fluxuations  in  Rb  concentration.  Initially, 
Rb  concentrations  are  below  optimal,  however  through  heating  the  vapor  pressure  is 
increased  until  the  most  effective  concentration  for  the  input  beam  strength  is  reached. 
Rollover  begins  as  alkali  concentration  becomes  high  and  the  gain  medium  optically 
thick.  Eventually  lasing  ceases  after  the  pump  beam  is  unable  to  penetrate  the  length  of 
the  gain  cell  and  gain  receeds  below  threshold.  This  design  shows  greater  control  over 
thermal  management  and  may  also  be  capable  of  being  expanded  to  a  time  varying 
temperature  control  system  that  is  based  on  the  current  output  efficiency  [17]. 

Since  their  advent  in  2003,  the  DEAL  has  seen  noteworthy  improvement  in  its 
realization.  Early  DEAL  demonstrations  achieved  maximum  output  powers  of  1  to  3 
Watts.  Currently  operating  CW  DEALs  are  in  the  lOO’s  of  Watts  and  the  pulsed  variety 
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in  the  lO’s  of  kilowatts.  Multiple  initiatives  into  solving  the  thermal  management  issues 
assoeiated  with  the  DPAL  have  shown  sueeess.  And  third,  impressive  slope  effieieneies 
have  been  observed  in  a  variety  of  experiments.  Research  and  design  is  ongoing  and  as 
more  data  is  collected  and  rates  and  quantities  defined  they  can  be  integrated  into  the 
model  making  it  more  precise  and  robust. 
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Ill,  Methodology 


Introduction 

Presented  here  is  a  derivation  from  fundamental  laser  physies  of  the  analytic  model 
developed  initially  by  Dr.  Gordon  Hager  and  extended  through  the  work  of  this  thesis. 
The  model  is  implemented  in  “Mathematica”,  a  leading  computation  software  package  in 
scientific  academia.  The  initial  derivation  will  neglect  pump  beam  linewidth  and  derive 
solutions  which  are  based  on  a  single  frequency  delta  function  input  beam.  Afterwards,  it 
will  be  shown  that  a  frequency  dependent  pump  beam  intensity  function  is  easily  applied 
to  the  solutions  of  this  more  simple  case.  The  model  incorporates  all  of  the  currently 
published  values  for  the  various  spectroscopic  constants  and  broadening  rates. 

While  this  model  is  analytic  in  nature,  it  does  rely  on  numerical  integration  of  certain 
functions,  and  some  approximations  are  made.  However,  these  functions  are  continuous 
over  the  bounds  of  integration  and  have  slopes  much  lower  than  the  numerical  integration 
step  size,  giving  confidence  in  resultant  values.  Also,  the  model  only  analyzes  a  Rb 
based  DPAL  with  an  ethane/helium  buffer  gas  mix.  Approximations  in  regards  to  certain 
effects  being  neglected  or  deemed  insubstantial  will  be  addressed  as  they  arise. 

As  mentioned  in  the  abstract,  the  model  produces  an  output  beam  intensity  function 
dependent  on  seven  input  values.  These  include  the  input  beam  intensity,  input  beam 
linewidth,  the  Rb  concentration,  ethane  concentration,  mirror  refiectivitiy,  cavity 
transmittivity,  and  gain  cell  length.  After  incorporating  all  the  fundamental  aspects  of  the 
physics  involved,  plots  are  generated  through  “Mathematica”  of  the  various  effects 
changes  in  system  parameters  have  on  output  beam  intensity.  The  figure  on  the 
following  page  highlights  all  the  essential  features  included  in  the  model. 
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T  ransmitti\'ity 


Figure  8.  Model  Dependent  Variables 


Shown  in  the  standard  layout  above,  the  model  has  yet  to  incorporate  spatial 
dynamics  and  is  independent  of  cavity  design.  It  is  currently  only  capable  of  predicting 
output  intensities,  not  net  output  power,  cavity  mode,  or  output  beam  quality.  Being  so, 
design  engineers  planning  a  system  based  on  the  results  of  the  model  will  have  freedom 
in  design,  constraining  themselves  only  to  meeting  the  parametric  settings  of  the  desired 
predicted  output  intensity.  The  pump  linewidth  extension  has  yet  to  be  made  on  the 
lasing  transition  and  its  cross-sectional  value  is  taken  at  its  linecenter  peak  value.  The 
table  on  the  following  page  is  to  be  used  as  reference  throughout  the  derivation  of  the 
model  and  lists  all  symbols  used  with  their  definition  and  value  (if  applicable). 
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Table  3.  Model  Variable  Definitions  (Rb  DPAL) 


Symbol 

Definition 

Value 

Units 

N, 

concentration  of  atoms  in  ground  state 

variable 

em'' 

N2 

eoneentration  of  atoms  in  the  5  P1/2  state 

variable 

em'^ 

N3 

eoneentration  of  atoms  in  the  5  P3/2  state 

variable 

em'' 

Ntot 

alkali  eoneentration 

variable 

em'' 

021 

absorption  eross  section  on  the  lasing  transition 

*4.80x10'*' 

(linecenter) 

em^ 

031 

absorption  eross  seetion  on  the  pump  transition 

*4.92x10'*' 

(lineeenter) 

em^ 

0 

longitudinally  averaged  two-way  intraeavity  pump  beam 

intensity 

variable 

W/cm^ 

h 

Planek’s  Constant 

6.626x10''^ 

J-s 

V21 

frequeney  of  the  laser  beam 

3.774x10*^ 

Hz 

V31 

frequeney  of  the  pump  beam 

3.846x10*^ 

Hz 

k32 

spin-orbit  mixing  rate  (Ethane  partner) 

5.2x10'*° 

cm  Vs 

t2I 

radiative  lifetime  of  upper  laser  level 

27.7 

ns 

t31 

radiative  lifetime  of  upper  pump  level 

26.23 

ns 

intraeavity  average  two-way  laser  intensity 

variable 

W/cm^ 

E 

Ethane  eoneentration 

variable 

cm'' 

4 

gain  eell  length 

variable 

cm 

R 

net  eavity  refleetivity 

variable 

N/A 

T 

net  eavity  transmission 

variable 

N/A 

K 

Boltzmann  eonstant 

1.38065x10'^^ 

J/K 

T 

Gain  eell  temperature 

variable 

K 

gth 

threshold  gain 

variable 

cm'* 

Ipin 

input  beam  intensity 

variable 

W/cm' 

c 

speed  of  light  (in  vacuum) 

3.0x10^ 

m/s 

hj 

pump  wavelength 

780.026 

nm 

hi 

laser  wavelength 

794.76 

nm 

Avsi 

pump  beam  bandwidth 

variable 

Hz 

He 

Helium  eoneentration 

variable 

cm'' 

gi 

degeneraey  of  ground  state 

2 

N/A 

g2 

degeneraey  of  upper  laser  level 

2 

N/A 

g3 

degeneraey  of  upper  pump  level 

4 

N/A 
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Longitudinally  Averaged  Two-Way  Intracavity  Pump  Beam  Intensity 

We  begin  by  defining  the  intracavity  pump  beam  intensity.  Consider  the  interaction 
between  the  pump  beam  at  each  transition  through,  or  reflection  off  of,  a  component  of 
the  cavity.  The  intensity  off  the  initial  interaction  with  the  beam  splitter  defines  Ipi„,  from 
there  the  beam  is  incident  on  the  first  wall  of  the  gain  cell,  receiving  a  transmission  loss 
defined  by  tp,  resulting  in  an  initial  inner  cell  intensity  of  tplpm.  As  the  beam  progresses 
through  the  gain  cell,  it  receives  a  Beer’s  Law  exponential  gain  given  by: 

(10) 

where  g3i  is  the  gain  on  the  pump  transition.  After  exiting  the  cell  another  transmission 
loss  is  taken  as  well  as  a  reflectivity  loss  defined  as  Vp.  Finally,  before  re-entering  the  cell 
it  experiences  another  transmission  loss  as  it  passes  through  the  cell  wall.  This  gives  the 
intensities  of  the  two  initial  input  beams  incoming  from  both  ends  of  the  cell  as: 

(11) 

(12) 

where  is  the  initial  incident  pump  beam  intensity  and  I~  is  the  intensity  of  the  input 
beam  after  a  single  pass  a  reflection  back  into  the  gain  medium.  Figure  9  below  shows 
the  input  beam  intensity  at  each  transition  for  a  single  round  trip  pass  [2]. 


Figure  9.  Round  Trip  Pump  Beam  Intensity 
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Assuming  the  beam  is  highly  attenuated  after  a  single  round  trip,  the  longitudinally 
averaged  intraeavity  pump  beam  intensity  ean  be  defined  by  integrating  the  sum  of  these 
two  ineident  beams  aeross  the  length  of  the  gain  medium,  and  averaging  the  result.  The 
explieit  derivation  is  shown  here; 

(13a) 

-  “  “  (13b) 

- - -  (13e) 

Substituting  values  from  equations  1 1  and  12  for  and  /“  yields: 

- - -  “  “  (13d) 

Ignoring  transmission  losses  and  assuming  100%  pump  beam  refleetivity  simplifies  this 
equation  to: 

- - -  “  (13e) 

Considering  the  beam  is  highly  attenuated  after  its  first  pass  and  the  first  transmission 
loss  eould  be  ineorporated  into  the  initial  beam  intensity,  the  approximation  in  eqn.  13e 
will  be  used  to  define  Q..  In  the  same  manner  output  laser  intensity  ean  be  defined  in 
terms  of  the  average  intraeavity  lasing  intensity  by  [2]; 


(14) 
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Laser  Rate  Equations 

The  dependence  of  both  the  intracavity  pump  beam  and  output  laser  intensity 
functions  on  alkali  population  densities  demand  solutions  to  the  dynamic  laser  rate 
equations  for  evaluation.  Assuming  no  second  order,  non-linear,  or  extraneous  effects  in 
the  gain  medium,  the  conservation  of  mass  and  the  rate  equations  for  the  three  states  in 
the  DPAL  during  steady  state  operation  are  defined  as  follows: 

—  -  -  -  -  (15) 

—  -  —  -  (16) 

—  -  “  -  (IV) 

(18) 

The  first  term  in  each  rate  equation,  and  the  second  term  in  eqn.  15  are  the  stimulated 
emission  and  absorption  rate  between  either  the  pump,  lasing  beam  or  both.  The  second 
term  in  eqns.  16  and  17  are  the  dynamics  of  the  spin-orbit  Boltzmann  equilibration. 
Finally,  the  remaining  terms  are  the  spontaneous  emissions. 

Solutions  to  these  equations  aren’t  readily  available  due  to  their  non-linearity, 
however,  eliminating  eqn.  17  still  gives  an  algebraic  set  which  is  easily  solved  through 
Mathematica.  There  are  two  sets  of  solutions,  one  for  the  non-lasing  mode  and  one  for 
lasing  mode.  In  the  non-lasing  mode  instance,  ¥  =  0  and  for  saturating  pump  intensities 
the  population  densities  are  defined  by  their  values  described  in  the  literature  review.  In 
lasing  mode,  'P  ^  0  and  since  the  gain  is  clamped  at  threshold  eqn  19  at  the  top  of  the 
following  page  can  be  used  to  solve  for  the  populations  and  intracavity  lasing  intensity 
simultaneously  [2]. 
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(19) 


The  simultaneous  solutions  to  the  the  steady  state  laser  mode  population  densities  and 
intraeavity  laser  beam  intensity  are  given  by: 

- — -  (20) 

- — -  (21) 

“  “  “  (22) 
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Transcendental  Solution 

We  now  have  solutions  to  the  independent  variables  (population  densities)  for  the 
intraeavity  pump  beam,  intraeavity  lasing  beam,  and  output  laser  intensity.  However, 
these  variables  are  themselves  dependent  on  the  averaged  intraeavity  pump  beam 
intensity.  This  eauses  the  solution  for  Q.  to  be  transeendental  in  nature,  as  in  it  is 
inherently  a  funetion  of  itself.  Rewriting  equation  13e  to  show  this  dependenee  explieitly: 

-  “  (24) 


27 


As  seen,  in  order  to  evaluate  and  obtain  values  for  the  intraeavity  pump  intensity  an 
input  value  for  said  intensity  must  determined.  It  is  only  when  both  the  independent 
variable  value  for  Q.  and  the  funetion  output  value  itself  are  equal  is  there  a  solution. 
These  are  found  by  solving  the  planar  intersection  between  f(f2)  =  and  The  graph 
below  is  an  example  of  this  intersection  for  typical  DPAL  parametric  settings. 


Figure  10.  Transcendental  Intersection 


The  intersection  of  the  two  functions  shown  above  are  the  solutions,  given  the 
baseline  parametric  settings,  to  the  intraeavity  longitudinally  average  pump  beam 
intensity  for  values  ranging  from  0  -  20000  W/cm  .  The  model  solves  this 
intersection  problem  by  solving  for  H  in  the  root  equation  below: 

-  “  -  (25) 
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In  doing  so  explicit  values  which  can  be  used  as  input  into  the  intracavity  laser 
intensity  equation  can  be  found  and  the  model  is  capable  of  predicting  output  laser  beam 
intensities  based  on  the  seven  fundamental  parameters  within  the  system.  Figure  1 1 
below  gives  Q.  values  for  a  range  of  input  beam  intensities  with  the  baseline  parametric 
settings  found  in  table  4. 


Figure  11.  £2vs/p„ 


Given  these  solutions  we  rewrite  the  output  beam  intensity  as  an  explicit  function  of 
input  beam  intensity,  and  thereby  completing  the  frequency  independent  part  of  the 
model.  With  this  ability,  we  now  have  a  “black  box”  capable  of  taking  in  a  set  of 
parametric  settings  and  giving  us  our  desired  result,  output  beam  intensity,  explicitly; 


(26) 
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Frequency  Dependence  of  the  Pump  Transition 

As  mentioned,  matching  the  pump  beam  spectral  profile  to  the  absorption  cross- 
section  lineshape  is  paramount  for  optimal  efficiency.  In  order  to  model  this,  we  must 
assign  a  frequency  dependence  to  the  absorption  cross-section  (031)  and  the  input  beam 
intensity  {Ipin).  This  is  done  by  appending  normalized  frequency  distribution  functions 
onto  the  mentioned  existing  values: 

(27) 

(28) 

Where  fpump  is  a  normalized  Gaussian  Distribution  and is  a  pressured  broadened 
Lorentzian  Distribution.  They  are  defined  explicitly  in  equations  29  and  30. 

- -  (29) 

-  -  (30) 

For  baseline  operating  conditions  of  600  Torr  of  Ethane  and  Rb  vapor  concentrations 


corresponding  to  393.15  K  the  pressure  broadened  width  of  the  absorption  cross  section 
has  a  FWHM  value  of  19.73  GHz.  By  setting  the  diode  FWHM  to  this  optimistically 
narrow  bandwidth  excellent  spectral  overlap  is  achieved  as  shown  in  figure  12  below: 


Fre^*encv  CH2> 

Figure  12.  D1  Spectral  Overlap  [3] 
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The  majority  of  the  produced  results  set  the  FWHM  of  the  absorption  cross  section  to 
a  predetermined  value  not  coupled  to  cell  temperature  or  buffer  gas  concentration. 
However  a  few  plots  incorporate  the  broadening  effects  associated  with  increased 
pressure  and  cell  temperature,  for  those  the  FWHM  is  defined  as: 

(31) 

Where  is  the  published  broadening  rate,  is  the  temperature  at  which  the  rate  was 
measured,  E  is  the  ethane  concentration  (or  corresponding  gas)  and  T  is  the  current 
DPAL  operating  temperature. 

Due  to  the  frequency  dependence  of  the  input  beam  and  absorption  cross-section,  the 
first  term  in  equation  15  becomes  an  integral  over  frequency  and  algaebraic  solutions  to 
the  rate  equations  cease.  However,  a  solution  to  the  problem  has  been  arrived  through 
redefining  the  rate  equations  with  an  overall  pump  beam  interaction  term  Qf.  It  begins  by 
adding  a  frequency  dependence  to  the  single  frequency  solution  of  D: 

- - -  “  (32) 

Multiplying  by  the  denominator  (except  for  /g)  and  integrating  over  frequency  gives: 


(33) 

The  term  on  the  left  of  eqn.  33  is  the  frequency  dependent  absorption  and  stimulated 
emission  rate  for  the  lasing  transition,  and  is  the  term  limiting  an  algebraic  solution.  By 
defining  Qf  as  this  and  substituting  it  into  the  rate  equations  solutions  are  once  again 
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obtained  for  the  population  densities,  however  this  time  they  are  dependent  on  Qf  This 
once  again  leads  to  a  transcendental  equation  for  Qf,  explicitly: 

(34) 

Solutions  are  again  obtained  through  the  root  equation  method,  and  new  population 
density  functions  are  generated  independent  of  Qf  and  are  now  functions  of  Ipi„.  By 
doing  so  the  “black  box”  is  once  again  established  and  the  frequency  addition  to  the 
pump  interactions  has  been  accomplished.  The  output  laser  intensity  is  now  a  function  of 
pump  beam  bandwidth  as  shown  in  eqn.  35. 

(35) 

With  the  ability  to  predict  output  laser  intensity  with  8  dimensions  of  variability  in  a 
quasi-analytic  fashion  the  DPAL  three  level  model  is  an  excellent  tool  for  analysis  of 
projected  future  technologic  endeavors.  The  few  approximations  that  have  been  made 
are  minimized  or  have  been  deemed  insubstantial.  For  a  complete  layout  of  the  code, 
with  a  brief  Mathematica  tutorial,  see  appendix  b. 
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IV,  Results  and  Analysis 


Introduction 

The  code  generates  2D  and  3D  plots  describing  laser  performance  characteristics 
including  efficiency,  pump  beam  output  intensity,  population  density  ratios,  direct 
comparisons  to  single  and  multi-frequency  cases,  and  others.  The  single  frequency 
results  are  presented  first  and  will  be  used  for  comparison  with  the  broadband  pump 
beam  results.  Each  plot  is  accompanied  with  an  analysis  of  the  results  given  the 
implications  that  the  graphs  trends  have  on  laser  characterization.  The  3D  plots  are 
capable  of  showing  the  simultaneous  effect  changes  in  multiple  systems  parameters  have 
on  the  laser  output  making  it  easier  to  assign  significance  to  each  input  variable. 

The  graphs  generated  are  based  on  a  set  of  “baseline”  parametric  settings.  These 
correspond  to  recent  experimental  operating  conditions  and  the  values  for  the  frequency 
independent  segment  of  the  code  are  given  in  the  table  below.  The  Rb  concentration 
corresponds  to  a  cell  temperature  of  393.15  K,  the  ethane  concentration  corresponds  to  a 
partial  pressure  of  300  Torr. 

Table  4.  Baseline  Parameters  -  Frequency  Independent  Case 


Variable 

Value 

Units 

h 

2 

cm 

^tot 

1.85x10^^ 

cm'^ 

E 

9.66x10^^ 

cm'^ 

r 

0.2 

N/A 

t 

0.975 

N/A 
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Single  Frequency,  No  Lasing  Results 

Consider  the  pre-lasing  conditions  of  the  DPAL  under  the  frequency  independent 
model.  The  steady  state  population  densities  as  a  function  of  fl  are  plotted  in  figure  13: 
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Figure  13.  Steady  State  Alkali  Concentrations  -  Pre  lasing 

As  seen,  as  intracavity  pump  intensity  increases  each  state  asymptotically 
approaches  their  values  described  above  in  the  literature  review. 

Figure  14  shows  the  intracavity  intensity  as  a  function  of  input  beam  strength. 
Initially  the  pump  strength  is  not  sufficient  to  fully  penetrate  the  gain  medium  and 
becomes  fully  absorbed  before  a  single  pass  through  the  cell.  This  can  be  seen  in  the  low 
input  beam  powers  of  the  graph  where  the  slope  efficiency  is  small.  However,  at 
approximately  225  W/cm  ,  the  beam  is  able  to  fully  pass  through  the  cell  and  a  portion 
of  the  reflected  wave  contributes  to  .  Eventually  the  gain  medium  becomes  bleached 
and  the  two  way  intracavity  beam  intensity  reaches  its  maximum  slope  efficiency  of  2. 

(36) 
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Single  Frequency,  Lasing  Mode  Results 

Presented  in  this  section  are  the  results  for  the  frequency  independent,  steady  state 
lasing  mode  case.  The  baseline  parameters  are  the  same  from  the  previous  section,  and  it 
is  shown  that  these  results  coincide  well  the  the  pre-lasing  mode  results.  Analysis  of  the 
frequency  independent  case  will  build  intuition  on  the  dynamics  and  interactions  of  the 
system  paving  the  way  for  a  clearer  view  of  the  effect  that  broadening  the  pump  beam  has 
on  the  DPAL.  Recall  in  this  instance  the  interaction  between  the  pump  beam  and  the 
system  is  represented  with  a  delta  function  in  frequency  input  intensity  centered  on  the 
peak  absorption  cross-section  value.  Being  so,  we  assume  this  represents  the  upper  limit 
on  absorption  and  for  certain  settings  efficiency,  output  beam  strength,  and  inversion 
ratios  for  a  given  input  beam  intensity.  Figure  15  on  the  following  page  shows  the 
oscillating  intracavity  laser  intensity  as  a  function  of  the  intracavity  beam  intensity. 
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Figure  15.  Intracavity  Laser  Intensity  -  Frequency  Independent 

Initially  the  pump  beam  is  unable  to  reach  threshold  and  'F  is  undefined.  However, 
at  an  Q.  value  of  40.04  W/cm  threshold  is  reached  and  lasing  begins.  The  linear  slope 
regime  above  represents  the  Q2L  behavior  with  the  spin-orbit  relaxation  able  to  keep  up 
with  the  pump  and  laser  transitions.  Rollover  occurs  as  these  radiative  transfer  rates 
dominate  the  dynamics  and  the  buffer  gas  is  unable  to  equilibrate  the  alkali  in  time  for  the 
relatively  high  oscillating  ¥  value  to  see  an  inversion.  An  asymptotic  limit  representing 
a  bleached  laser  transition  begins  to  develop  at  higher  Q.  values.  This  threshold  value 
corresponds  well  with  the  calculated  threshold  value  of  40.3  W/cm  for  the  no  lasing 
case,  showing  the  model  is  well  coordinated  and  precise. 
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The  plot  of  the  intraeavity  pump  beam  intensity  as  a  function  of  input  beam  strength 
for  the  lasing  case  is  similar  to  the  non  lasing  except  the  curvature  is  described  by 
different  phenomena.  As  seen  below,  once  again  H  is  slow  to  initially  rise,  then  turns  and 
hits  a  maximum  slope  efficiency.  Initially  the  spin-orbit  rate  is  able  of  cycling  the  atoms 
into  the  upper  lasing  state  where  the  inversion  with  ground  maintains  amplification  of  the 
oscillating  laser  intensity.  However,  for  the  defined  parametric  settings,  absorption 
begins  to  cease  as  the  spin-orbit  rate  is  unable  to  funnel  atoms  into  the  lower  fine- 
structure  level.  This  results  in  a  linear  slope  efficiency  of  over  100%.  This  plot  is  shown 
in  Figure  16  below: 


Figure  16.  Intracavity  Pump  Intensity  -  Frequency  Independent 
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The  dramatic  effect  that  varying  ethane  concentration  has  on  output  intensity  and 


slope  efficiency  is  readily  seen  in  Figure  17  shown  here: 


Figure  17.  Ethane  Variation  -  Frequency  Independent 
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As  seen,  once  the  maximum  spin-orbit  relaxation  rate  is  reached  rapid  rollover 
begins  as  the  cyclic  photon  engine  begins  to  break  down,  and  bottlenecking  in  the  upper 
pump  state  occurs.  While  increasing  the  buffer  gas  pressure  seems  like  the  apparent 
solution,  broadening  rates  on  the  cross-sections  must  be  considered  and  maintaining  a 
single  mode  laser  intensity  becomes  more  challenging  as  more  and  more  modes  rise 
above  threshold.  There  is  also  potential  for  cross-section  overlap  at  high  enough 
pressures  which  would  greatly  reduce  efficiency  of  the  ability  to  lase. 

As  an  inverse,  consider  holding  the  ethane  pressure  a  constant  and  varying  alkali 
concentration.  In  doing  so  similar  results  are  produced  as  can  be  seen  on  Figure  18 
presented  here: 
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Figure  18.  Rubidium  Variatiou  -  Frequeucy  ludepeudeut 


When  the  alkali  number  density  is  low,  the  eell  beeomes  easily  bleaehed  by  the 
ineident  pump  beam  and  the  majority  of  the  ineident  photons  are  wasted.  However,  the 
Q2L  limit  is  maintained  at  higher  input  beam  intensities  and  a  linear  slope  effieiency 
maintained.  The  eoneentrations  are  not  so  optieally  thiek  as  to  limit  lasing  at  sueh  high 
pump  intensities. 

Window  transmittivity  also  plays  a  vital  role  in  slope  effieieney.  Gain  is  dependent 
on  the  negative  of  the  log  of  the  transmission  factor,  implying  that  as  transmittivity 
increases,  gain  becomes  more  negative. 
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Figure  19.  Transmittivity  Variation  -  Frequency  Independent 


It  is  evident  that  Brewster  windows  and  other  transmission  boosting  design  elements 
will  play  a  eritieal  role  in  optieal-optieal  effieieney. 

Cell  Length  is  easily  eharaeterized  and  the  optimal  length  is  where  the  final  round 
trip  value  of approaehes  zero.  Figure  20  at  the  top  of  the  following  page  shows  the 
rapid  rollover  in  slope  effieieney  as  the  eell  beeomes  bleaehed  by  high  pump  intensities. 

It  shows  that  the  longer  the  eell,  the  longer  the  slope  effieieney  stays  linear  and  relatively 
high.  This  is  attributed  to  the  idea  that  the  pump  beam  beeomes  more  dispersed 
throughout  a  larger  medium  and  therefore  less  radiation  is  ineident  per  photon.  This 
reduees  absorption  on  the  pump  transition,  allowing  for  the  buffer  gas  to  stay  ahead  of  the 
game  and  within  the  Q2L  regime. 
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Figure  20.  Cell  Length  Variation  -  Frequency  Independent 


Broadband  Pump,  Lasing  Mode  Results 

With  confidence  in  the  single  frequency  model  the  addition  of  pump  bandwidth  has 
been  made  and  the  results  are  shown  here.  Consider  first  only  varying  pump  beam 
bandwidth.  It  is  apparent  from  Figure  21  on  the  following  page  that  matching  the 
bandwidth  to  the  absorption  has  a  strong  effect  on  optical-optical  efficiency.  It  also  has 
effects  on  the  rollover  point,  population  inversion  and  3D  plots  showing  the  effects  of 
bandwidth  on  the  other  parameters  are  analyzed.  Comparison  to  the  single  frequency 
case  shows  a  reduction  in  efficiency  as  expected.  The  graph  on  the  following  page  gives 
our  frequency  dependent  pump  beam  intensity  as  a  function  oilpm  for  a  range  of  spin- 
orbit  pressures. 
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Figure  21.  Q  vs  Ip/„  -  Frequency  Dependent 
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The  highest  bandwidth  (500  GHz),  has  relatively  low  slope  efficiency  as  a  somewhat 
small  proportion  of  the  photons  are  within  the  absorption  profile  of  the  transition.  As 
bandwidth  decreases,  slope  efficiency  rises.  This  is  caused  by  the  increase  in  absorption 
and  a  more  rapid  approach  to  the  bleached  transition  limit.  As  you  can  see,  the  10  GHz 
pump  beam  mimics  the  single  frequency  case  in  that  around  5200  W/cm  the  cell 
becomes  bleached  and  approaches  its  limiting  value  of  21  pm.  However,  all  other 
bandwidths  excluding  the  500  Ghz  have  intracavity  pump  beam  intensities  greater  than 
Ipin,  implying  that  absorption  towards  a  saturated  state  is  consistent  with  those 
bandwidths.  Given  that  more  broadened  pump  beams  can  saturate  the  D2  transition,  there 
may  be  less  of  a  demand  to  pursue  narrowbanded  diode  stacks  .  An  investigation  into 
overall  efficiency  between  narrowbanding  the  diode  bars  or  pumping  broadband  may 
prove  worthwhile.  A  second  plot  emphasizing  the  importance  of  spectral  overlap  is  seen 
at  the  top  of  the  following  page. 
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Figure  22.  I;„,^  vs  Ip„,  -  Frequency  Dependent 

Figure  22  illustrates  the  importance  of  spectral  overlap.  This  chart  was  produced 
with  600  Torr  of  ethane  in  order  to  ensure  the  results  stayed  in  the  Q2L  limit.  It  shows 
the  output  laser  intensity  as  a  function  of  input  beam  intensity  with  FWHM’s  ranging 
from  the  absorption  cross-sections  FWHM  of  10  GHz,  up  to  200  GHz.  As  is  apparent, 
without  a  tight  Gaussian  pump  beam  efficiency  is  drastically  reduced  and  net  output 
power  greatly  diminishes.  As  bandwidth  increases,  more  and  more  photons  lie  outside  of 
a  significant  absorption  probability,  causing  not  only  wasted  incident  radiation  but 
polluting  the  cavity  with  extraneous  energy  potentially  leading  to  thermal  aberrations  in 
the  gain  medium. 

As  bandwidth  increases  and  absorption  declines  the  Q2L  limit  is  extended  beyond  its 
single  frequency  point.  This  is  most  evident  through  a  comparison  of  Figures  22  and  23 
on  the  following  page  and  Figure  17.  They  are  the  same  plots  except  Figure  17  is  the 
single  frequency  case.  Figure  22  has  a  Gaussian  pump  beam  with  a  FWHM  of  19.73 
GHz,  while  Figure  23  has  a  FWHM  of  200  GHz. 
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Figure  23.  Ethane  Variation  -  Frequency  Dependent  (19.73  GHz) 


Shown  above  are  the  effeets  of  ethane  eoneentration  on  output  intensity  for  a  19.73 
GHz  input  beam.  Through  eomparison  with  Figure  17,  it  is  shown  that  the  linear  slope 
effieieney  is  maintained  for  larger  values  of  I  pin.  This  is  due  to  the  reduetion  in  absorbed 
beam  intensity  as  bandwidth  is  inereased.  This  can  be  seen  to  an  even  greater  extent  in 
Figure  24  below  (200  GHz  input  beam  bandwidth): 


Figure  24.  Ethane  Variation  -  Frequency  Dependent  (200  GHz) 
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The  reflectivity  of  the  cavity  at  the  lasing  wavelength  seems  to  play  a  rather 
insignificant  role  in  output  beam  intensity.  This  may  be  attributed  the  extremely  high 
gain  nature  of  the  alkali  gain  medium  and  its  ability  to  reach  threshold  regardless  of 
cavity  reflection.  The  only  obvious  yet  significant  effect  is  seen  when  r  =  1.0,  or 
complete  reflection  back  into  the  cavity.  This  result  strengthens  the  confidence  in  the 
model  as  another  intuitive  result  is  produced  as  expected.  Figure  25  show  output 
intensities  for  reflectivities  ranging  from  0.2  -  0.5.  The  zero  output  intensity  straddles 
the  axis  and  is  unseen  but  present  on  a  color  graph. 
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Figure  25.  Reflectivity  Variation  -  Frequency  Dependent 


As  seen,  reflectivity  within  the  cavity  slightly  effects  slope  efficiency.  It  does  not  seem 
to  effect  the  rollover  effect  or  even  maximize  output  beam  intensity.  This  strengthens  the 
idea  that  the  key  rate  or  dynamic  of  the  system  is  the  spin-orbit  relaxation  rate. 
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With  this  in  mind  a  new  parameter,  A,  is  defined  as  N3/N2.  It  represents  how  elose 
the  upper  two  laser  levels  are  to  their  Boltzmann  Distribution,  or  how  close  the  system  is 
to  operating  in  the  Q2L  regime.  To  see  how  the  pump  bandwidth  effects  this  ratio,  a  3D 
contour  plot  with  iSspump  and  Ipin  as  the  independent  variables  is  generated  and  shown 
below.  The  upper  bound  to  the  Boltzmann  Equilibrated  ratio  for  baseline  parametric 
settings  is  set  at  A  =  1.73. 


Figure  26.  A  -Bandwidth  Variation  -  Frequency  Dependent 


As  seen  pump  bandwidth  plays  a  significant  role  in  defining  A.  At  lower 
bandwidths,  A  rises  more  quickly  as  expected  through  a  higher  net  absorption.  It 
continues  to  stay  above  higher  bandwidth  values  for  the  same  pump  intensity  but  rises 
more  slowly  as  A  approaches  its  upper  limit.  This  may  be  attributed  to  the  decreased 
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percentage  of  the  populations  out  of  equilibrium  which  may  reduce  the  buffer  gas 
effective  equilibrium  rate.  At  high  pump  bandwidths  the  rise  in  A  is  linear  corresponding 
to  the  Q2L  regime.  In  this  case  as  more  photons  are  absorbed  into  the  upper  lasing  level 
there  are  more  opportunities  for  the  buffer  gas  to  equilibrate  with  the  alkalis  in  this  state. 
It  also  shows  that  this  rate  is  able  to  keep  up  with  the  stimulated  emission  rate  on  the 
lasing  transition. 

While  the  2D  plot  of  reflectivity  suggested  a  weak  dependence  between  output 
intensity  and  reflectivity,  the  contour  plot  below  suggests  that  a  lower  reflectivity 
produces  a  higher  output  laser  intensity.  It  shows  that  as  reflectivity  drops  output 
intensity  rises  in  a  non-linear  fashion. 
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Figure  27.  Reflectivity  Variation  3D  -  Frequency  Dependent 


47 


As  seen  maximum  laser  intensity  is  aehieved  at  zero  refleetivity.  This  onee  again 
suggest  alkalis  possess  extremely  high  gain  and  a  single  ASE  output  beam  maybe  be 
eapable  of  produeing  higher  intensities  than  an  oscillating  laser  cavity.  However,  for 
high  beam  quality  a  proper  single  mode  resonator  should  be  designed  with  a  relatively 
low  reflectivity. 

Parametric  efficiency  studies  have  been  conducted  generating  DPAL  efficiency  plots 
for  a  range  of  input  beam  strengths,  buffer  gas  concentrations  and  pump  beam  bandwidth. 
Figure  25  below  represents  the  optical-optical  efficiency  for  a  DPAL  operating  with  a  10 
kW  input  beam,  300  Torr  of  Ethane  and  50  Torr  He.  The  results  incorporate  the 
broadening  rates  recently  published  by  Pitz  et  al. 


DPAL  BGciency 


Figure  28.  DPAL  Efficiency  -  Frequency  Dependent 


DPAL  efficiency  is  plotted  against  the  natural  log  of  the  Rb  concentration  and  the 
ratio  between  the  pump  and  absorption  cross  section  bandwidth.  As  expected,  efficiency 
depends  highly  on  bandwidth.  This  graph  also  shows  the  effects  of  high  Rb  concentration 
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on  output.  At  the  lower  Rb  eoneentrations,  the  pump  beam  is  fully  absorbed  and  slope 
effieiencies  are  linear  and  high  for  eaeh  given  bandwidth.  Near  the  center  value  of  the 
Rb  concentration  efficiency  drops  as  the  gain  medium  becomes  optically  thick  for  a  10 
kW  input  beam.  This  rollover  however  occurs  at  higher  and  higher  Rb  concentrations  as 
bandwidth  in  increased,  as  expected. 

As  mentioned  above  reflectivity  plays  a  key  role  in  determining  output  beam 
strength.  Presented  below  are  the  output  intensities  for  a  80  kW  input  beam  at  various 
bandwidths.  Once  again,  as  bandwidth  rises,  output  drops. 
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This  plot  reinforces  the  results  of  the  3D  reflectivitiy  contour  plot  above  in  that 
output  intensity  rises  as  reflectivity  drops.  However,  eventually  at  high  enough  pump 
intensities  these  curves  converge  into  a  single  line.  This  probably  represents  a  pump 
beam  intensity  that  bleaches  the  pump  transition  at  all  bandwidths.  Figure  30  on  the 
following  page  gives  the  output  intensity  vs  reflectivity  for  a  100  kW  input  beam,  and 
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shows  this  effect.  The  model  uses  a  baseline  value  for  cavity  reflectivity  of  0.2, 
corresponding  to  an  efficiently  designed  laser  system.  At  r  =  0.2,  there  is  still  high  output 
intensities  but  with  the  ability  to  potentially  create  a  resonating,  single  mode,  diffraction 
limited  output  beam. 
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As  seen,  at  high  pump  intensities,  the  effects  of  bandwidth  become  insignificant  due  to 
saturation.  This  effectively  reduces  the  system  to  is  frequency  independent  model  and 
solutions  converge  upon  the  lowest  bandwidth  instance. 

One  particular  value  of  interest  when  analyzing  pump  bandwidth  effects  on  laser 
performance  is  the  ratio  between  the  frequency  independent  intracavity  pump  beam 
intensity  (O)  and  the  frequency  dependent  intracavity  pump  beam  intensity  (Of). 
Explicitly,  Gamma  is  defined  as; 

-  (36) 
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Intuitively  one  would  believe  that  the  frequeney  independent  intraeavity  pump  beam 
intensity  is  bounded  to  be  always  greater  than  the  frequency  dependent  case.  However, 
upon  execution  of  the  code  the  results  are  somewhat  counter-intuitive.  Figure  3 1  Gamma 
at  a  variety  of  pump  bandwidths  for  the  baseline  parametric  settings. 
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And  the  same  plot  again  for  a  more  pressurized  (500  Torr)  gain  cell: 
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Figure  32.  Gamma  -  500  Torr  Ethane 
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Figure  33  -  O  vs  Ipi„  500  Torr  Ethane  (500  Ghz) 


Figure  33  shows  the  eause  for  Gamma  dropping  below  1.  Initially  the  broadband  pump 
has  a  greater  slope  effieiency,  caused  by  a  reduction  in  absorption.  However,  at  higher 
pump  intensities  and  saturated  transition  conditions  the  single  frequency  case  approaches 
is  intracavity  slope  limit  of  2Ip/„,  while  the  broadband  case  maintains  a  lower  slope  .  At 
around  15000  W/cm  the  single  frequency  rate  has  bleached  the  transition  wheeras  the 
broadband  case  has  not,  leading  to  a  greater  slope  efficiency.  Returning  again  to  laser 
intensity.  Figure  34  shows  output  intensity  as  a  function  of  Rb  concentration  for  a  set  of 
pump  bandwidths: 


Figure  34.  IiaseVS.  Log[Rb]  -  500  Torr  Ethaue 
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As  is  apparent  in  the  graph,  there  is  an  optimal  eombination  of  alkali  eoneentration, 
and  bandwidth  whieh  maximizes  output  and  this  eombination  points  direetly  at  low 
bandwidth.  At  any  given  Rb  eoneentration,  the  lowest  bandwidth  has  the  highest  output 
intensity,  and  at  10  GHz  the  appropriate  Rb  eoneentration  results  in  an  optieal  to  optieal 
effieieney  of  about  80%.  A  seeond  plot  of  F  vs  y  for  a  family  of  pump  bandwidths  is 
shown  at  the  top  of  the  following  page.  In  order  for  the  pump  beam  to  drive  the  system  it 
must  exeite  quieker  than  the  lasing  beam  relaxes.  In  order  to  do  so  it  must  get  absorbed 
whieh  is  evident  in  Figure  34.  At  lower  bandwidths  A  maintains  an  inversion  and  there  is 
bottleneeking  at  the  upper  pump  state.  This  is  a  potential  souree  of  inefficieney  and 
further  studies  of  a  A  for  max  effieieney  may  prove  benefieial. 


Figure  35.  A  -  300  Torr  Ethane 
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As  shown,  A  deereases  as  the  spin  orbit  rate  inereases,  or  as  Ethane  is  added  to  the 
system.  A  low  A  is  good  for  effieieney  but  one  that  is  too  low  may  limit  laser 
performanee. 
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The  model  is  eapable  of  produeing  a  variety  of  results  that  can  be  used  to  analyze 
and  more  intuitively  understand  the  dynamics  and  various  interactions  of  the  system.  By 
accurately  predicting  the  output  beam  intensity  based  on  the  fundamental  elements  of  any 
laser  system,  parametric  solutions  can  be  obtained  upon  which  engineers  can  base  their 
designs. 
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V,  Conclusions  and  Discussion 


With  the  results  of  the  model  readily  available  and  interpretable,  conclusions  on  the 
effects  of  broadband  pumping  are  made.  Also  discussed  are  various  effects  between  the 
other  independent  variables.  Third,  issues  with  the  approximations  made  are  addressed 
and  finally,  suggestions  for  future  work  on  the  model  are  made. 

In  essence,  a  broadening  of  the  pump  bandwidth  offers  very  little  towards  achieving 
maximum  laser  performance.  As  seen,  as  bandwidth  increases,  output  laser  intensity 
decreases.  Given  that  an  objective  of  the  DPAL  is  cost  efficiency  and  low  power-usage, 
maximizing  pump  absorption  is  critical.  However,  the  model  does  show  that  a  very  high 
efficiency  is  theoretically  possible,  (albeit  an  engineering  feat  yet  to  be  seen).  There  are 
also  certain  cases  where  it  may  be  possible  for  a  more  broad  banded  beam  to  actually  be 
more  efficient,  as  is  seen  in  Figure  21. 

While  it  may  seem  intuitive,  in  a  quantum  mechanical  sense,  that  any  broadening  of 
an  energizing  spectral  source  will  reduce  efficiency,  the  current  state  of  affairs  regarding 
diode  stacks  puts  their  bandwidths  in  the  lOO’s  of  GHz.  With  this  in  mind,  if  demand  for 
output  power  begins  to  heavily  outweigh  efficiency  concerns,  systems  can  be  designed 
based  on  the  appropriate  pump  beam  interaction,  allowing  the  model  to  “pave  the  way” 
for  future  progress.  The  model  is  also  an  open  book  without  an  ending  and  is  capable  of 
become  far  more  robust,  accurate,  precise  and  efficient. 

As  seen  throughout  the  data,  there  must  be  a  merger  between  all  system  parameters 
in  order  to  maximize  performance.  However,  it  was  shown  that  some  variables  have  a 
greater  influence  than  others.  For  example,  reflectivity  plays  a  minor  role  in  determining 
overall  performance  whereas  buffer  gas  concentrations  are  far  more  critical.  Bandwidth 
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plays  a  significant  role  however  in  praetieality  it  seems  as  if  it  would  be  easier  to  simply 
inerease  the  pump  intensity  rather  than  adjust  the  ethane  eoneentration.  Cell  length  in  the 
standard  DPAL  design  is  held  eonstant  and  the  various  other  parameters  ean  be  adjusted 
to  make  the  most  out  of  any  length  eell.  In  this  author’s  opinion,  the  system  primarily 
eomes  down  to  alkali  and  buffer  gas  eoneentrations,  and  pump  beam  intensity. 

While  the  model  is  quasi-analytie,  founded  in  laser  engineering,  and  eompares  well 
with  experimental  data,  it  lacks  a  sense  of  eompleteness  or  of  an  all  eneompassing 
simulation  of  a  true  laser  system.  In  a  way,  it  is  an  algorithm  redueed  to  the  key  elements 
of  any  laser  system,  a  simplifieation,  an  “ideal”  representation  of  a  far  more  detailed 
phenomenon.  Of  greatest  eoneem  is  the  idea  that  all  the  alkalis  reside  in  one  of  only 
three  possible  states  (  for  Rb:  5  P3/2,  5  P1/2  or  5  S1/2).  Reeently  eollected  data  suggests 
that  alkali’s,  upon  being  exeited  to  the  5  P3/2  state  experienee  a  myriad  of  kinetie 
proeesses  and  exeitations  before  fluoreseing  back  to  ground.  This  results  in  “spreading 
out”  of  the  alkalis  atoms  into  the  various  eleetron  states.  Fluoreseence  of  a  laser  pumped, 
vacuum  sealed,  thermally  insulated  Rb  eell  was  eollected  by  means  of  an  optical  cable 
and  a  0.3  M  speetrometer.  A  series  of  data  eolleetions  were  undertaken  to  observe  the 
intensity  fiuetuations  of  the  various  Rb  emission  lines  as  a  funetion  of  Rb  eoneentration 
and  ineident  beam  intensity.  For  eaeh  set  the  Rb  vapor  was  subjeet  to  high  intensity, 
pulsed  radiation  in  resonanee  with  the  D2  transition.  A  Grotrain  Diagram  displaying  the 
various  observed  emission  lines  is  shown  on  the  following  page. 
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Figure  36.  Grotrain  Diagram  -  Laser  Pumped  Rb  Fluorescence 
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As  one  can  see,  many  different  energy  levels  become  populated  solely  from  pumping 
on  the  D2  line.  A  closer  inspection  of  Figure  29  reveals  that  all  dipole  allowed  transitions 
are  observed  and  even  some  non-dipole  allowed.  A  detector  with  a  larger  bandwidth  may 
reveal  even  more  Rb  emissions. 

A  second  issue  with  the  model  is  that  it  makes  great  strides  in  an  effort  to  maintain  its 
quasi-analytic  nature.  However,  at  times,  this  author  feels  too  much  of  the  physical 
phenomena  occurring  in  the  laser  system  is  ignored  or  deemed  insignificant.  In  the 
derivations  of  the  longitudinally  averaged  two  way  intracavity  pump  beam  intensity 
certain  terms  are  ignored  or  neglected  for  simplicity.  This  author  suggests  that  one 
should  take  advantage  of  the  computational  power  accessible  by  AFIT  and  other  facilities 
to  produce  a  more  thorough  account  of  the  interactions  within  the  system. 

Even  so,  the  model  still  makes  accurate  predictions  and  encompasses  all  of  the 
fundamental  laser  physics  involved  with  the  DPAL.  It  is  an  excellent  tool  for  analysis, 
exploration,  brainstorming,  and  there  are  many  upgrades  that  could  begin  tightening  the 
gap  begin  the  model  and  the  lab.  First,  extending  the  frequency  dependence  of  the  pump 
transition  to  the  lasing  transition  would  produce  more  accurate  and  fundamentally  sound 
results.  Second,  considering  numerical  integration  is  already  implemented  in  the  code,  a 
more  precise  and  approximation  free  definition  of  the  various  system  variables  may  open 
up  new  phenomena  or  unforeseen  issues  within  the  system.  For  user  -friendliness, 
implementing  the  code  into  a  compiled  language  would  dramatically  reduce  runtime. 

The  closer  the  model  represents  reality,  the  closer  it  will  predict  it. 
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Appendix  A 


Rubidium  Physical  Data: 


Electronic  data 

Shells:  2,8,18,8,1 
Orbitals:  [Kr]  5s^ 
Electronegativity:  0.9,  0.8 

1.  Ionization  potential: 

4.1771  eV 

2.  Ionization  potential: 

27.28  eV 

3.  Ionization  potential:  40.0 
eV 

Oxidation  states:  1 


Physical  data 

Thermal  data 

Melting  point:  38.89  °C 
Boiling  point:  686  °C 
Specific  heat:  0.363  J/gK 
Heat  of  fusion:  2.192  kJ/mol 
Heat  of  vaporization:  72.216 
kJ/mol 

Thermal  conductivity:  0.582 
W/cmK 


Steric  data 

Atomic  radius:  2.98  A 
Ionic  radius:  1.61  A 
Covalent  radius:  2.16  A 
Atomic  volume:  55.9  cm^/mol 
Density  (293  K):  1 .53  g/cm^ 
Crystal  structure:  Cubic:  Body 
centered 


Isotopes 


Nuclide 

Abundance  [%] 

Mass 

Spin 

Halflife 

Decay  mode 

®'Rb 

0 

81 

3/2 

4.57h 

EC,/S^ 

®^Rb 

0 

82 

1 

1.258m 

P\EC 

“Rb 

0 

83 

5/2 

86.2d 

EC,EC 

«^Rb 

0 

83.914 

2 

32.9d 

EC,P\P~ 

®®Rb 

72.15 

84.9118 

5/2 

-- 

Stable 

“'’Rbm 

0 

85 

0 

>1ps 

IT 

®®Rb 

0 

85.911 

2 

18.65d 

P~,EC 

®^Rb 

27.85 

86.9092 

3/2 

4.8E10y 

P- 

®®Rb 

0 

88 

2 

17.7m 

P- 

®®Rb 

0 

89 

3/2 

15.44m 

P- 

®°Rb 

0 

90 

1 

2.6m 

P- 

®°Rbm 

0 

90 

4 

4.3m 

P',\T 
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J^pendix  B  -  Simulation  Code  (Three  Level  Laser  System) 


Presented  here  is  a  Brief  Mathematica  Tutorial  and  a  concisely  annotated 
computer  simulated  model  for  The  Diode  Pumped  Alkali  Laser. 

Mathematica  Tutorial  : 

This  tutorial  will  focus  on  the  necessities  of  executing  the  code . 

1 .  )  Cell  Evaluation  : 

On  the  far  right  of  each  window  are  small  "Cell"  dividers  and  allow 
one  to  select  a  cell  by  left  clicking  on  the  divider.  Once  selected, 

hit  either  the  enter  key  on  the  number  pad  or  shift  + 
enter.  Do  not  evaluate  cells  containing  notation. 

2 .  )  Execution  of  Code  : 

a.  )  First  select  a  section  of  code  tailored  to  your  analysis  goals 
(i.e.  -  if  you  want  to  do  frequency  independent  analysis, 

select  one  of  the  frequency  independent  sections  of  the  appendix  below) 

b .  )  Once  selected ,  input  all  user  dependent  variables  into  the 
"System  Parameters" subsection. 

c.  )  After  con^letely  defining  all  parametric  values, 

evaluate  all  cells  not  consisting  of  notation. 

d .  )  At  the  end  of  each  section  is  a  rubrick  giving  an  outline  of  how  to 
generate  plots  corresponding  to  that  section.  Singly  input  the  user 
specified  plot  range  and  independent  varickble  and  execute  the  cell. 

3 .  )  Generating  Black  and  White  Plots  with  Markers  and  Legends  : 

To  generate  sophisticated  plots  capable  of  describing  more  data 
simultaneously  utilize  the  "MultipleListPlot"  function.  This  function 
is  capable  of  generating  marker  points  for  each  data  point  euid  a  legend 
f or  each  data  set . 

a.)  Load  graphics  package  (ignore error  message)  ; 

«  Graphics 'MultipleListPlot ' ; 
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b. )  Rubrick  for  generating  Black  and  White  Plots  of  multiple  lines  : 


MultipleListPlot[ 

{Tablel  [{X  -  varicible,  Y  -  Variable}  ,  {X,  X  -  min,  X  -  max,  step -size}], 

Table2  [{X  -  variable,  Y  -  Variable} ,  {X,  X  -  min,  X  -  max,  step -size}]. 

Tables  [{X  -  variable,  Y  -  Variable]  ,  {X,  X  -  min,  X  -  rneix,  step-size}], 

Table4  [{X  -  variable,  Y  -  Variable} ,  {X,  X  -  min,  X  -  max,  step  -  size}] } , 

PlotLabel  -►  "Insert  Title",  FrameLabel  -►  {"X-Axis  Label",  "Y-Axis  Label"}, 
GridLines  -»  Automatic,  Frame  True,  Fillings tyle  -»  None, 

Plot  Joined  -►  True, 

PlotLegend  -►  {Tablel  label,  Table2  label.  Tables  label,  Table4  label}] 

4.)  Exporting  data  as  a  .dat  file  (ASCII) 

To  generate  a  file  consisting  of  a  table  of  the  independent  and  dependent 
variables  of  a  function  utilize  the  "Export"  function. 

Exan^le  : 

i .  )  Generate  a  Table  : 

Tablel  =  Table  [{Ipin,  15*Ipin},  {Ipin,  0,  500,  10}] 

{{0,  0],  {10,  150],  {20,  SOO],  {SO,  450],  {40,  600],  {50,  750],  {60,  900], 

{70,  1050},  {80,  1200},  {90,  1350},  {100,  1500},  {110,  1650},  {120,  1800}, 

{ISO,  1950},  {140,  2100},  {150,  2250},  {160,  2400},  {170,  2550},  {180,  2700}, 

{190,  2850},  {200,  3000},  {210,  3150},  {220,  3300},  {230,  3450},  {240,  3600}, 

{250,  3750},  {260,  3900},  {270,  4050},  {280,  4200},  {290,  4350},  {300,  4500}, 

{310,  4650},  {320,  4800},  {330,  4950},  {340,  5100},  {350,  5250}, 

{360,  5400},  {370,  5550},  {380,  5700},  {390,  5850},  {400,  6000}, 

{410,  6150},  {420,  6300},  {430,  6450},  {440,  6600},  {450,  6750}, 

{460,  6900},  {470,  7050},  {480,  7200},  {490,  7350},  {500,  7500}} 

ii.  )  Export  the  tcible  specifying  save  location  and  file  type  : 

Export {"C: WDocuments  cuid  SettingsWpjonesWTablel.dat",  Tablel] 
C:\Document3  and  Setting3\pjone3\Tablel.dat 

The  file  can  be  fovind  at  the  file  path  given  by  the  output. 

This  concludes  the  tutorial . 
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**  NOTE  *♦ 


1 .  Each  section  is  to  be  viewed  as  cm  independently  operating  simulation.  Within 
each  part  is  a  con^lete  system  model . 

2 .  Upon  execution  of  any  section  ensure  that  all  variables  have  been  cleared 
by  quitting  the  local  Mathematica  Kernel .  This  is  accon^lished  by 
selecting  the  "Evaluation"  menu  button  in  Mathematica  and  then  selecting 
"Quit  Kernel . " 

3 .  User  input  variables  have  been  cleared  and  inputs  are  to  be  specified 
under  "System  Parameters." 


Ethane  Number  Densities  (u/cm^)  Corresponding  to  Various  Pressures  (Torr)  : 


100  Torr 

3.22  *10“ 

200  Torr 

6.44*10“ 

300  Torr 

9.66*10“ 

400  Torr 

1.28*10^’ 

500  Torr 

1.61*10“ 
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Baseline  Parametric  Settings  : 


(Singly  execute  this  segment  if  you  wish  to  evaluate  the  code  at  its  baseline 
operating  conditions ,  emd  spectroscopic  constants  corresponding  the 
the  table  found  in  Reference  [] ) 

h  :=  6.626068*10'^^ 

t21  :=  27.70*10'® 

t31  :=  26.23*10'® 

021  :=  4.80*10'^® 

o31  :=  4.92*10'^® 

k32  ;=  3.16*10^-10 

v31  :=  3.845*10''14 

v21  :=  3.774*10^14 

e  :=  0.8670379156330899 

gth[r_,  t_,  igr_]  :=  -Log[r  *  t®]  /  (2 *  ig) 

k  :=  1.3806504  *10''-23 

T  :=  300 

c  :=  3.0*10'‘10 

A31  :  =  c  /  v31 

A21  :  =  c  /  v21 

M  :=  9.66*10^18 

ntot  :=  1.85*10''13 

Ig  :=  2 

r  :=  0.2 

t  :=  0.975 

r  :=  M*k32 
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Parti:  Frequency  Independent,  No  Lasing 


(Assume  steady  state) 

Rate  Equations  ($  =  0) 

Solve  f 

n3 

t31 ' 


0  ==  nl  +  n2  +  n3  -  ntot, 
0  ==  $}, 


r  Q  $  n2 

1 0  ==  031  *  (n3  -  2  *  nl)  * -  +  a21  *  (n2  -  nl)  * -  +  -  + 


h»v31 


h*v21  t21 


S  n2 

0  ==  -021  *  (n2  -  nl)  * -  +  If  *  (n3  -  2  *Exp[-e]  *n2) - , 

h*v21  t21 


{nl,  n2,  n3}j 

This  solves  for  our  population  densities  as  a  function  of  intracavity  pun^ 
beam  intensity  : 

nl[i?J  :  = 

(ntot  (e®  h  v31  +  2  h  t21  v31  y  +  e®  h  t31  v31  y  +  e®  t31  o31  J2  +  2  t21 131  y  o31  S) )  / 
(e®  h  v31  +  2  h  t21  v31  y  +  e®  h  t31  v31  y  +  3  c®  t31  o31  +  6  t21 131  y  o31  J2  + 

2  e®  t21 131  y  o31  J2) 

n2[J2_]  :=  (2  c®  ntot  t21 131  y  o31  J2)  / 

(e®  h  v31  +  2  h  t21  v31  y  +  c®  h  t31  v31  y  +  3  c®  t31  o31  i?  +  6  t21 131  y  o31  5  + 

2  e®  t21 131  y  o31  O) 

n3[i^]  :=  (2  ntot  t31  (e®  +  2  t21  y)  o31  i2)  / 

(e®  h  v31  +  2  h  t21  v31  y  +  c®  h  t31  v31  y  +  3  c®  t31  o31  J2  +  6  t21 131  y  o31  i?  + 

2  e®  t21 131  y  o31  O) 
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System  Parameters  : 
h  :=  6.626068*10'^* 
t21  :  = 
t31  :  = 

021  :  = 

031  :  = 
k32  :  = 
v31  :  = 
v21  :  = 

d  :=  0.8670379156330899 

gth[r^,  t_,  igr_]  :=  -Log[r  *  t*]  /  (2  *  icr) 

k  :=  1.3806504  *10'‘-23 

T  :  = 

c  :=  3.0*10^10 
A31  :=  c/v31 
A21  :=  c/v21 
M  :  = 
ntot  : = 

Ig  :  = 
r  :  = 

t  :  = 

r  :=  M*k32 
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Defining  our  Two  -  Way  Intracavity  Longitudinally  Averaged  Pun^Beam  Intensity 
(Frequency  Independent) 

i'  Ipin 

QQ[Ipin  ,  Q  ]  :=  -  * 

“  “  la31*  (n3[P]  -  2*nl[P])  *lg.' 

(Exp [031*  (n3[i?]  -  2*nl[i?])  *2*lg]  -  1) 

Q[Ipin_]  :=  (0  /.  FindRoot[«fi[IpijQ,  Q]  -  Q,  {Q,  (Ipin)}]) 

Rubrick  for  generating  plots  : 

Plot[Q[Ipin] ,  {Ipin,  Minimum,  Maximum},  GridLines  -»  Automatic, 

Frame  -*  True,  PlotLabel  "Insert  Title  Here", 

FrameLabel  -»  {"X-Axis  Label"  ,  "Y-Axis  Label"}] 
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Part  2  :  Frequency  Independent ,  Lasing 
(Assume  steady  state) 

Rate  Equations  (S  #  0) 


Solve 


|o  ==  031*  (n3  -  2  *nl)  * 


$  n2  n3 


+  021  *  (n2  -  nl)  * 


h*v31 


n2 


0  XX  -o21*  (n2  -  nl)  * 


+  y*  (n3  -  2*Exp[-©]  *n2) 


h*v21 


t21 


0  x=  nl  +  n2  +  n3  -  ntot, 

-Log[r*t*]  /  (2*lg)  =x  o21  *  (n2  -  nl)J, 

{nl,  n2,  n3,  S}J  //  FullSimplify 

Defining  our  steady  state  lasing  populations  as  functions  of  our  input 
parameters  : 

nl[i^,  iitot_,  M_,  r_,  t_,  icr_]  :  = 

(2hk32Mt31v31Log[r  t*]  +  e®  (h  (v31  +  k32  M  t31  v31)  +  t31  o31  J2) 

(2  La  ntot  o21  +  Log[r  t^] ) )  / 

(4  ig  o21  (h  k32  M  t31  v31  +  c®  (h  (v31  +  k32  M  t31  v31)  +  2  t31  o31  S) ) ) 

n2[j;^,  ntot_,  M_,  r_,  t_,  icr_]  :  = 

(c®  (2  ig  ntot  o21  (h  (v31  +  k32  M  t31  v31)  +  t31  o31  O)  - 
(h  (v31  +  k32  M  t31  v31)  +  3  t31  o31  J2)  Log[r  t*]))  / 

(4  ig  o21  (h  k32  M  t31  v31  +  e®  (h  (v31  +  k32  M  t31  v31)  +  2  t31  o31  S) ) ) 

n3[i^,  ntot_,  M_,  r_,  t_,  ig_]  :  = 

(t31  (2  Lantota21  (h k32  M  v31  +  e®  o31  J2)  +  (-h k32  M  v31  +  e®  o31  j2)  Log[r  t®]))  / 
(2  ig  o21  (h  k32  M  t31  v31  +  c®  (h  (v31  +  k32  M  t31  v31)  +  2  t31  o31  O) ) ) 
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Defining  our  Two  -  Way  Intracavity  Longitudinally  Averaged  Pur^Beam  Intensity  : 
(Frequency  Independent) 

Q[Ipin_,  ntot_,  M_,  r_,  t_,  lgr_]  :  = 

(Q  /.  FindRoot[$iQ[Ipin,  Q,  ntot,  M,  r,  t,  Ig]  -  Q,  {Q,  (Ipin)}]) 

QQ[Ipin_,  S_,  ntot_,  M_,  r_,  t_,  icr_]  :  = 

(Ipin  /  (031*  ((n3[i3,  ntot,  M,  r,  t,  ig])  -  2*  (nl[g,  ntot,  M,  r,  t,  ig]))  •  ig))  * 
(Exp[a31  *  ( (n3  [S,  ntot,  M,  r,  t,  ig] )  -  2*  (nl  [j?,  ntot,  M,  r,  t,  ig] ) )  *  2  *  ig]  - 
1) 

Our  Intracavity  Laser  Intensity  : 

$[Ipin_,  ntot_,  M_,  r_,  t_,  ig_]  :  = 

-  (h  v21 

(-2  ig  ntot 021  (2  k32  M  t21  (hv31  +  t31o31  *  (Q[Ipin,  ntot,  M,  r,  t,  ig]))  + 

(h  (v31  +  k32  M  t31  v31)  +  (1  -  2  k32  M  t21)  t31  o31  * 

(Q[Ipin,  ntot,  M,  r,  t,  ig])))  + 

(2  k32  M  t21  (h  v31  +  3  t31  o31*  (Q[Ipin,  ntot,  M,  r,  t,  ig]))  + 
e®  (h  (v31  +  k32  M  t31  v31)  +  (3  +  2  k32  M  t21)  t31  o31  * 

(Q[Ipin,  ntot,  M,  r,  t,  ig] )) )  Log[r  t^] ) )  / 

(2  t21  o21 

(hk32  M  t31  v31  + 

€®  (h  (v31  +  k32  M  t31  v31)  +2t31o31*  (0[lpin,  ntot,  M,  r,  t,  ig]))) 

Log[r  t^]) 


Next  we  Define  our  Output  Laser  Intensity  : 

Ilias*  ntot_,  H_,  r_,  t_,  ig_]  :  = 

$[Ipin,  ntot,  M,  r,  t,  ig]  * 

(gth[r,  t,  ig]  *  ig *Exp[gth[r,  t,  ig]  *  ig]  *  t*  (1  -  r))  / 

( (Exp[gth[r,  t,  ig]  »ig]-l)*(l  +  t^*  r»Exp[gth[r,  t,  ig]  *  ig] ) ) 
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Rubricks  for  generating  plots  : 


□iDega  Plots  : 

Plot[Q[lpin,  1.85*10^^  9.66*10'‘18,  0.2,  0.975,  2] 
GridLines  -»  Automatic,  Frame  True,  PlotLabel 
FrameLcibel  -►  {"X-Axis  LcQjel"  ,  "Y-Axis  Label"}] 

Laser  Plots  : 

Plot[{IIi»„[Ipin,  ntot,  3.22*10^18,  r,  t,  Ig]  , 
IIiaie[Ipin,  ntot,  6.44*10'‘18,  r,  t,  Ig]  , 
IIiaseEIpin,  ntot,  9.66*10‘'18,  r,  t,  Ig]  , 
IIiaie[Ipin»  ntot,  1.28*10'^19,  r,  t,  Ig]  , 
IIias«[Ipin,  ntot,  1.61*10‘'19,  r,  t,  Ig]},  {Ipin, 
GridLines  -»  Automatic,  Frame  -»  True,  PlotLabel  -» 
FrameLabel  -+  {"X-Axis  Label"  ,  "Y-Axis  Label"}] 


,  {Ipin,  0,  12000}, 
"Insert  Title  Here", 


0,  12  000}, 

"Insert  Title  Here", 
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Parts:  Frequency  Dependent,  Lasing 
(Assume  steady  state) 

Rate  Equations  (S  *  0) 


Solve  ^ 
{O  == 


-Qf  5  n2  nS 

-  +  021  *  (n2  -  nl)  * -  +  -  +  - , 

h*v31  h*v21  t21  tSl 


0 


$ 

-021*  (n2  -  nl)  *  -  +  k32*M*  (n3  -  2*Exp[-0]  *n2) 

h*v21 


n2 

t21 ' 


0  ==  nl  +  n2  +  n3  -  ntot, 


-Log[r*t*]  /  (2*lg)  ==  o21  *  (n2  -  nl)|, 

{nl,  n2,  n3,  5}]  //  FullSin^lify 

Defining  our  state  population  densities  and  intracavity  laser  intensity  : 

nl[i2^,  ntot_,  M_,  r_,  t_,  icr_]  :  = 

(2  e®  ig  o21  (hntot  (1  +  k32  M  t31)  v31  -  t31  Qf)  + 
h  (2  k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  Log[r  t*] )  / 

(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31)  )  v31  o2l) 

n2  [i2f_,  iitot_,  M_,  r_,  t_,  ig_]  :  = 

(e®  (2  ig  o21  (hntot  (1  +  k32  M  t31)  v31  -  t31£(f)  -  h  (1  +  k32  M  t31)  v31  Log[r  t^]))  / 
(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  o2l) 

n3[j2f_,  ntot_,  M_,  r_,  t_,  ig_]  :  = 

(2hk32  ig  M  ntot  t31  v31  o21  +  2  c®  ig  t31  o21  J2f  -  h  k32  M  t31  v31  Log[r  t®] )  / 

(2  €®  h  ig  v31  o21  +  2  h  k32  ig  M  t31  v31  o21  +  2  e®  h  k32  ig  M  t31  v31  o2l) 

5[Ipin_,  ntot_,  M_,  He_,  r_,  t_,  ig_,  4vpuii?j_]  :  = 

-  (v21 

(4  k32  ig  M  t21  o21 

(-hntot  v31  +  t31  *  (QQf[Ipin,  ntot,  M,  He,  r,  t,  ig,  ^vpun^j]))  - 
2  e®  ig  o21  (h  ntot  (1  +  k32  M  t31)  v31  - 

(l  +  2k32Mt21)  t31*  (S)S}f[Ipin,  ntot,  M,  He,  r,  t,  ig,  ^vpuii?)]))  + 
h  (2k32Mt21  +  €®  (l  +  k32Mt31))  v31  Log[r  t®] ) )  / 

(2  t21  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  o21  Log[r  t^] ) 
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System  Parameters  : 
h  :=  6.626068*10'^* 
t21  :  = 
t31  :  = 

021  :  = 

031  :  = 
k32  :  = 
v31  :  = 
v21  :  = 

d  :=  0.8670379156330899 

gth[r^,  t_,  igr_]  :=  -Log[r  *  t*]  /  (2  *  icr) 

k  :=  1.3806504  *10'‘-23 

T  :  = 

c  :=  3.0*10^10 
A31  :=  c/v31 
A21  :=  c/v21 
M  :  = 
ntot  : = 

Ig  :  = 
r  :  = 

t  :  = 

r  :=  M*k32 
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Pressure  Dependence  of  Absorption  cross  -  section  (courtesy  Greg  Pitz) 

Av31[M_,  He_]  := 

(((((M)  *15*7*  (10^6)  *  (7.5006»10*-3))  »26.2)  *10^6)  + 

((((He)  *k*T*  (10*6)  *  (7 . 5006  *  10*-3) )  *20.0)  *10*6))  *  ( (T  /  310)  *  (1  /  2)  )  + 
5.2*10*6 


Cross  -  section  and  piu^  beam  lineshapes  : 


Avpun?j_] 


■2*Log[2]^/2  ' 
*  Avpimp  I 


*Exp[-4  *  Log[2]  *  ( (v  -  v31)  /Avpun?))^] 


0031  [V  ,  M  ,  He  ] 


(A3i)^ 

-  *  f31[v,  M,  He] 

8  *  TT  *  t31  I 

I 


f31[v  ,  M  ,  He  ] 


'  Av31[M,  He] 

i  2  *  TT  1 

((V  -  v31)^  +  (Av31[M,  He]  /2)^) 

Next  We  define  our  Input  Beam  Interaction  with  the  system  : 

QSJQf  [Ipin_,  i2f_,  ntot_,  M_,  He_,  r_,  t_,  icr_,  Avpiin5>_]  :  = 

N'lpin  *  fp.™  [v,  Avpuii5>]  ' 

- I  * 

(1- 

Exp[ 

(0031 [V,  M,  He] * 

(n3[i2f,  ntot,  M,  r,  t,  ig]  -  2*nl[i2f,  ntot,  M,  r,  t,  ig] )  *  2  *  ig)  ] ) , 
{v,  v31  -  10^^  v31  +  10“-®}] 

QQf[Ipin_,  iitot_,  H_,  He_,  r_,  t_,  ig_,  Avpuii?>_]  :  = 

(Qf  /.  FindRoot[QQQf [Ipin,  Qf,  ntot,  M,  He,  r,  t,  ig,  AvpulI^?]  -  Qf, 

{Qf,  (Ipin)}]) 
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Redefining  our  population  densities  as  a  function  of  Ipin  : 

nnl[Ipin_,  ntot_,  M_,  He_,  r_,  t_,  ig_,  4vpuH?)_]  :  = 

(2  e®  Ig  021 

(h  ntot  (1  +  k32  M  t31)  v31  - 

t31  *  ($}l}f[Ipin,  ntot,  M,  He,  r,  t,  ig,  4vpuii?)]))  + 
h  (2  k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  Log[r  t*] )  / 

(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  a2l) 

nn2[Ipin_,  iitot_,  H_,  He_,  r_,  t_,  ig_,  ,dvpuii5)_]  :  = 

(2  ig  o21  (h  ntot  (1  +  k32  M  t31)  v31  - 

t31  *  (Qflf[Ipin,  ntot,  M,  He,  r,  t,  ig,  ^ivpuiQ?]))  - 
h  (1  +  k32  M  t31)  v31  Log[r  t*] ) )  / 

(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  a2l) 

nn3[Ipin_,  ntot_,  He_,  r_,  t_,  ig_,  ,lvpini?j_]  :  = 

(2  h  k32  ig  M  ntot  t31  v31  o21  + 

2  e®  ig  t31  a21  *  (QQf  [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpu]i3>] )  - 
hk32Mt31v31Log[rt*])  / 

(2  e®  h  ig  v31  o21  +  2  h  k32  ig  M  t31  v31  a21  +  2  e®  h  k32  ig  M  t31  v31  a2l) 

Finally  our  Laser  Output  as  a  function  of  our  system  parameters  : 

IIlase[Ipin_,  ntot_,  M_,  He_,  r_,  t_,  ig_,  nvpini®i_]  :  = 

(a21  *  {nn2[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun^?]  - 

nnl [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun?)] )*ig*t*(l  -r)* 
Exp[o21*  (nn2[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun^?]  - 
nnl  [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpump])  *  ig]  )  j 
((Exp[a21*  (nn2[Ipin,  ntot,  M,  He,  r,  t,  la,  nvpia^]  - 

nnl  [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpuii?)] )  »  ig]  -  1)  » 

(1  - 

t^  *  r  * 

Exp[a21*  (nn2[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun^?]  - 
nnl  [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpunp] )  *  ig] ) )  * 
®[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpui®)] 
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Rxibrick  for  Generating  Plots  : 


In  this  case  it  is  suggested  to  create  plots  based  on  lists  of  data  due  to 
the  high  computational  demands  of  evaluating  the  frequency  dependent 
output  beam  intensity  : 

MultipleListPlot [ 

{Table[{lpin,  IIlase[lpin,  1.85*10“,  9.66*10^°,  1.61*10^18,  0.2, 

0.975,  2,  10*10''9]},  {Ipin,  0,  25000,  750}], 

Table[{lpin,  IIlase[lpin,  1.85*10^^,  9.66*10“,  1.61*10^18,  0.2, 

0.975,2,  50*10^9]},  {Ipin,  0,  25000,  750}], 

Table[{lpin,  IIlase[lpin,  1.85*10^^,  9.66*10“,  1.61*10^18,  0.2, 

0.975,  2,  100*10''9]},  {Ipin,  0,  25000,  750}], 

Table  [{Ipin,  IIlase[lpin,  1.85*10^^,  9.66*10“,  1.61  *10 '>18,  0.2, 

0.975,2,  200*10^9]},  {Ipin,  0,  25000,750}], 

Table[{lpin,  IIlase[lpin,  1.85*10“,  9.66*10^',  1.61*10'‘18,  0.2, 

0.975,  2,  500*10''9]],  {Ipin,  0,  25000,  750}]],  Gridlines  Automatic, 
Frame  -*  True,  PlotLabel  -*  "  lias*  vs  Ipir.  (Frequency  Dependent,  Lasing)", 
Fillings tyle  -»  None,  Plot Joined  -*  True, 

FrameLabel  ("Ipir  Watts/cm^",  "Ii„.  Watts/cm^"}, 

PlotLegend  -»  {lOGHz,  50GHz,  lOOGHz,  200  GHz,  500  GHz}] 
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Part  4  -  Gamma  Plots 


Presented  below  is  a  con®)lete  algorithmic  derivation  for  Gamma .  In  order 
to  generate  Gamma  plots  one  must  singly  input  the  user  pareuueters  and 
execute  the  code  below  : 

Defining  Omega  : 

Single  Frequency  : 

nl[J^,  ntot_,  M_,  He_,  r_,  t_,  ig_]  :  = 

(2ht31  v31  *  (k32»M)  *  Log[r  t^]  + 

e®  (h  (v31  +  t31  v31  *  (k32  *  M) )  +  t31  *  (o31)  *  S?)  (2  Lg  ntot  a21  +  Log[r  t*] ) )  / 

(4  Ig  021  (h  t31  v31  •  (k32  *  M)  +  c®  (h  (v31  +  t31  v31  *  (k32  *  M) )  +  2  t31  *  (a31)  *  1?) ) ) 

n2[i^,  ntot_,  M_,  He_,  r_,  t_,  ig_]  :  = 

(e®  (2  ig  ntot  a21  (h  (v31  +  t31  v31  *  (k32  *  M) )  +  t31  *  (a31)  *  i?)  - 
(h  (v31  +  t31  v31  *  (k32  *  M) )  +  3  t31  *  (a31)  *  i?)  Log[r  t*] ) )  / 

(4  ig  021  (h  t31  v31  *  (k32  «  M)  +  e®  (h  (v31  +  t31  v31  *  (k32  *  M) )  +  2  t31  *  (o31)  *  i?) ) ) 

n3[i^,  ntot_,  M_,  He_,  r_,  t_,  ig_]  :  = 

(t31  (2  igntoto21  (hv31*  (k32*M)  +  e®  *  (o31)  *i?)  + 

(-h  v31  *  (k32  *  M)  +  «®  *  (o31)  *  i?)  Log[r  t^] ) )  / 

(2  ig  o21 

(h  t31  v31  *  (k32  *  M)  +  e®  (h  (v31  +  t31  v31  *  (k32  *  M) )  +  2  t31  *  (o31)  *  i?) ) ) 

QS}fi[Ipin_,  I?_,  ntot_,  M_,  He_,  r_,  t_,  ig_]  ;  = 

(Ipin/  ((o31)  *  (n3ti?,  ntot,  M,  He,  r,  t,  ig]  -  2*nl[P,  ntot,  M,  He,  r,  t,  ig]  )  » 
ia))  * 

(Exp[(o31)  »  (n3[P,  ntot,  M,  He,  r,  t,  ig]  -  2*nl[j?,  ntot,  M,  He,  r,  t,  ig]  )  • 
2*  ig]  -  1) 

fll[Ipin_,  ntot_,  M_,  He_,  r_,  t_,  ig_]  :  = 

(Qf  /.  FindRoot[JM}Q[Ipin,  Qf,  ntot,  M,  He,  r,  t,  ig]  -  Qf,  {Qf,  (Ipin)}]) 


75 


Broadband  Case  : 


nbl[J2^,  ntot_,  M_,  r_,  t_,  igr_]  :  = 

(2  e®  in  021  (hntot  (1  +  k32  M  t31)  v31  -  t31  i2f)  + 
h  (2  k32  M  t31  +  c®  (1  +  k32  M  t31) )  v31  Log[r  t*] )  / 

(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  a2l) 

nnbl[Ipin_,  ntot_,  He_,  r_,  t_,  Lg_,  4vpiii5)_]  :  = 

(2  e®  ig  o21 

(h  ntot  (1  +  k32  M  t31)  v31  - 

t31  *  (QQf[Ipin,  ntot ,  M,  He,  r,  t,  ig,  idvpini^?]))  + 
h  (2  k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  Log[r  t®] )  / 

(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  a2l) 

nb2[J2^,  ntot_,  M_,  r_,  t_,  ig_]  :  = 

(e®  (2  ig  o21  (hntot  (1  +  k32  M  t31)  v31  -  t31i2f)  -h  (1  +  k32  M  t31)  v31  Log[r  t^]))  / 
(4  h  ig  (k32  M  t31  +  e®  (1  +  k32  M  t31) )  v31  a2l) 

nnb2[Ipin_,  ntot_,  1^,  He_,  r_,  t_,  ig_,  4vpiii5)_]  :  = 

(2  ig  o21  (h  ntot  (1  +  k32  M  t31)  v31  - 

t31»  ($M2f[Ipin,  ntot,  M,  He,  r,  t,  ig,  ^vpun?>]))  - 
h  (1  +  k32  M  t31)  v31  Log[r  t*]))  / 

(4  h  ig  (k32  M  t31  +  c®  (1  +  k32  M  t31) )  v31  a2l) 

nb3[j2f_,  ntot_,  H_,  r_,  t_,  ig_]  :  = 

(2hk32  ig  M  ntot  t31  v31  o21  +  2  e®  ig  t31  o21  J2f  -  h  k32  M  t31  v31  Log[r  t*] )  / 

(2  e®  h  ig  v31  a21  +  2  h  k32  ig  M  t31  v31  a21  +  2  e®  h  k32  ig  M  t31  v31  o2l) 

nnb3[Ipin_,  ntot_,  M_,  He_,  r_,  t_,  ig_,  ,lvpai5)_]  :  = 

(2  h  k32  ig  M  ntot  t31  v31  o21  + 

2  e®  ig  t31  021  *  (QQf[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpuz^?])  - 
hk32  M  t31  v31  Log[r  t*])  / 

(2  e®  h  ig  v31  o21  +  2  h  k32  ig  M  t31  v31  o21  +  2  e®  h  k32  ig  M  t31  v31  o2l) 
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QQQf  [Ipin_,  Of_,  Btot_,  M_,  He_,  r_,  t_,  igr_,  4vpini?)_]  :  = 

N'lpin  *  [v,  Jvpun?)]  ' 

- I  * 

(1- 

Exp[ 

(aa31[v,  M,  He]  * 

(nb3[fif,  Btot,  M,  r,  t.  La]  -  2*nbl[i2f,  utot,  M,  r,  t,  ig] )  *2  *  ig) ] ) , 
{v,  v31  -  v31  + 

Ql}f[Ipin_,  iitot_,  H_,  He_,  r_,  t_,  ig_,  4vpun?)_]  :  = 

(Qf  /.  FindRoot[QQQf [Ipin,  Qf,  utot,  M,  He,  r,  t,  Ig,  4vpuii?>]  -  Qf, 

{Qf,  (Ipin)}]) 

QQQQf[Ipin_,  iitot_,  M_,  He_,  r_,  t_,  ig_,  4vpuiig3_]  :  = 

NIntegrate [ 

((Ipin*  fpai?)[v,  ^vpump])  / 

(aCT31[v,  M,  He]  *  (rmb3[IplQ,  ntot,  M,  He,  r,  t,  Ig,  ^vpisq?]  - 
2  *nnbl  [Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun??]  )  *  ig) )  * 

(Exp[aa31[v,  M,  He] * 

(nnb3[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpump]  - 
2  » nnbl[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpungj]  )  »2  *  ig]  -  1) , 

(v,  v31  -  10^,  v31  +  10^}] 

Defining  Gamma  : 

r[Ipin_,  ntot_,  He_,  r_,  t_,  ig_,  nvpuii5)_]  :  = 

QQQQf[Ipin,  ntot,  M,  He,  r,  t,  ig,  nvpun®)]  /  Ql[Ipin,  ntot,  M,  He,  r,  t,  ig] 
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System  Parameters  : 
h  :=  6.626068*10'^* 
t21  :  = 
t31  :  = 

021  :  = 

031  :  = 
k32  :  = 
v31  :  = 
v21  :  = 

d  :=  0.8670379156330899 

gth[r^,  t_,  igr_]  :=  -Log[r  *  t*]  /  (2  *  icr) 

k  :=  1.3806504  *10'‘-23 

T  :  = 

c  :=  3.0*10^10 
A31  :=  c/v31 
A21  :=  c/v21 
M  :  = 
ntot  : = 

Ig  :  = 
r  :  = 

t  :  = 

r  :=  M*k32 
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Pressure  Dependence  of  Absorption  cross  -  section  (courtesy  Greg  Pitz) 

Av31[M_,  He_]  := 

(((((M)  *k»T*  (10^6)  *  (7.5006*10^-3))  *26.2)  *10*6)  + 

((((He)  *k*T*  (10*6)  *  (7 . 5006  *  10*-3) )  *20.0)  *10*6))  *  (  (T  /  310)  *  (1  /  2)  )  + 
5.2*10*6 


Cross  -  section  and  pun^  beam  lineshapes  : 

'2*Log[2]^/2 


fp.jin>[v  ,  Avpun?)  ]  :  = 


0031  [V  ,  M  ,  He  ]  :  = 


f31[v  ,  M  ,  He  ]  := 


*  Avpump 
'  (A3i)2  ' 


*Exp[-4  *  Log[2]  *  ( (v  -  v31)  /Avpun?))^] 


8  *  7r  *  t31 


*f31[v,  M,  He] 


1  '  1* 

'  Av31[M,  He] 

i  2  *  TT  1 

((v-v31)^+  (Av31[M,  He] /2)^) 
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